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•LOtlAIV  or  TCINI 


•nimmtk  - the  directlcm  of  afproech  of  the  eouad  e«»e  front,  deteralneS  by  the 
Crqjeetioo  of  the  ««re  froet  aorael  on  the  horiiootal  plene  tangent 
to  tiM  earth  at  the  point  of  observation. 

iiffrmction  - deviation  of  sound  raves  froe  a noraal  course  in  an  unobstructed 
aediua  when  paeaing  the  angle  or  edge  of  a large  body. 

dri/t-*scr«a«st  - the  total  dlsplaceeent  of  the  sound  rave,  noraal  to  its 

direction  of  propagation,  within  the  layer  of  ataosphere  considered. 

dri/t<4»tsd  - the  horiiootal  wind  coaponent  noraal  to  the  direction  of  sound 
propagation. 

dri/ttwtsd  frodisst  - the  rate  of  change  with  elevation  of  horiiontal  wind  coa- 
patent  noraal  to  the  direction  of  sound  propagation. 

iseiissttos  the  angle  of  the  sound  wave  noraal  with  the  horizontal  plane 
tangent  to  the  earth  at  the  point  of  observation. 

low  of  rtfrmeiion  - the  relation  connecting  the  velocity  of  sound,  wind,  and 
inclination  of  the  sound  rave  noraal  at  one  point  with  their  values 
at  another. 

loess  - a nuaber  of  possible  positions  satisfying  a given  set  of  conditions. 

Paraastsr  - a quantity  fixed  for  a given  set  of  aeteorological  conditions,  but 
varying  froa  one  set  of  conditions  to  another. 

Pormmttrie  solution  - a solution  in  teras  of  the  variables  and  a paraaeter, 
the  paraaeter  being  constant  once  the  set  of  conditions  is  defined. 

rsays'tlscrsasst  - the  total  displaceaent  of  the  sound  rave,  in  its  direction 
of  propagation,  within  the  layer  of  ataosphere  considered. 

rs/raetl*s  wind  - the  horiiontal  wind  coaponent  in  the  direction  of  sound 
propagation. 
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rgfrmetiv  ttind  gradUnt  - the  rate  of  change  eith  elevation  of  the  horlsontal 
wind  conponent  In  the  direction  of  sound  propagation. 


sound  path  - space  curve  traversed  by  a point  on  a sound  wave  surface. 
sound  ray  - 


sound  vsloeity  gradisnt  - the  rate  of  change  of  the  velocity  of  sound  with 
elevation. 


s^aes  array  - any  srrangenent  of  nicrophones  such  that  all  are  not  in  the 
sane  plane. 


sPteifie  humidity  - the  mass  of  eater  vapor  in  unit  nass  of  noist  air 


tstrahsdron  - a solid  bounded  by  four  plane  triangular  faces. 

tiiM-iscrvevst  - the  time  of  travel  of  the  sound  wave  along  a particular  path 
through  the  atnospheric  layer  considered. 

virtwai  tsmpsraturs  - the  temperature  of  absolutely  dry  air  having  the  sane 
density  as  moist  air  at  a slightly  lower  temperature. 
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Btratificstlon  parMotor,  defined  by  ratio  of  sound  velocity 
gradient  to  refractive  wind  gradient 

direction  angles  of  normal  to  sound  vave  front 

ratio  of  specific  heat  at  constant  pressure  to  specific  heat 
at  constant  volume 

Y as  defined  above  for  dry  air,  humid  (or  moist)  air,  and 
rater  vapor  respectively 

third  order  determinant 

a small  positive  quantity;  angular  correction 

true  azimuth  angle  from  north  (positive  clockwise) 

cosine  of  angle  of  Inclination  of  wave  front  normal  to  horizontal 
plane 

azimuth  of  sound  approach 

direction  cosines  of  normal  to  sound  wave  front 
180  degrees 
density 

density  of  dry  air 
density  of  raUr  vapor 
drift  psrsaeter 
angle  between  two  lines 
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lower  Md  upper  boundaries,  respectively,  of  s layer  of 
ataoephere 

velocity  of  sound  in  a still  ateosphere 
constant 

specific  heat  at  constant  pressure 
specific  heat  at  constant  voluae 
third  order  deteminant 

index  of  reliability  (shortest  distance  between  two  lines) 

total  differentials  with  respect  to  tine  and  space  coordinates, 
respectively 

sound  velocity  gradient 

virtual  tenperature  gradient  (or  lapse  rate) 

refractive  wind  gradient 

drift-wind  gradient 

partial  differentiation 

infinitesinal  increnents  in  space  coordinates 

tiae-increnent 

range-inereaent 

drlft-increasnt 

depth  of  ataospheric  layer 

e 


rinM&i 
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k coBstMt  of  integration;  range  paraaeter 

■4  aolecttlar  weight  of  dry  air 

■1,  aolecttlar  weight  of  water  vapor 


f point  in  space;  distance  in  s-direction  between  two  coordinate 

systeas 

P,M,i  coordinates  of  origin  of  one  coordinate  system  in  relation  to 

another 

p distance  of  plane  from  origin  of  coordinate  system;  atmospheric 

pressure 

P4  ataospherid  pressure  of  dry  air 

q specific  huaidity 

saturation  specific  humidity  (at  given  T and  p) 

i gas  constant,  equal  to  8.8703  X 10*  c.g.s  units  for  dry  air; 

distance  between  sound  source  and  observation  point 

i'  universal  gas  constant 

I 

r relative  huaidity 

Vg  distance  between  two  points 

f tiae  of  arrival  of  sound  wave;  temperature  of  air  in  degrees 

absolute 

V virtual  teaperature  in  degrees  absolute 
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vind  coaponents  In  x,  y,  x directions  respectively 
refractive  wind 

instantaneous  velocity  of  propagation  of  sound,  motion  of 
medium  included 

apparent  velocity  of  sound  along  the  horizontal,  constant  for 
a particular  sound  ray 

drift-wind 

grid-coordinates  of  point  of  observation 
grid-coordinates  of  sound  source 
integral  of  dx^  and  dy^.  respectively 
space  coordinates 

coordinates  after  clockwise  rotation  about  Z-axis,  of  amount 

- of  previous  coordinate  system 
2 

level  of  total  refraction 
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k naory  of  SowM  fropiptloii  ThroMgk  tia  M—pfcrt  aa  an  Applicatlai  to  Soaid  RMpiof 


I.  INTIOIVCTION 

This  eoMtltutes  Voluae  I of  s series  of  reports  dealing  vith  (a)  the 
effect  of  aeteorologlcal  conditions  on  the  propagation  of  sound  through  the 
ataosphere,  and  (b)  the  development  of  a technique  of  sound  ranging  whereby 
existent  meteorological  conditions  can  be  quantitatively  evaluated  in  terms 
of  their  net  effect  on  sound  propagation  from  source  to  microphone  position. 
Specifically,  it  is  the  intent  of  the  present  report  to  formulate  a theory  of 
sound  propagation  from  basic  analytical  considerations,  adequate  for  any  meteor- 
ological situation  whatsoever,  and  so  presented  that  it  nay  readily  be  applied 
to  military  sound  ranging  on  arti-llery. " 

An  adequate,  though  not  necessarily  feasible,  method  of  approach  to 
the  problem  of  the  determination  of  meteorological  corrections  for  sound  rang- 
ing would  be  one  concerned  with  an  analysis  of  sound  propagation  aiong  the 


particular  path  traversed  by  the  sound  wave.  The  concept  of  the  sound  path, 
or  ray,  is  not  new.  for  example,  Hilne^  had  occasion  to  investigate  an  analogous 
problem,  that  of  the  location  of  aircraft  by  means  of  sound  propagation  through 
the  atmosphere.  Gutenberg*  and  others  have  been  concerned  with  the  abnormal 
propagation  of  sound  to  great  distances,  and  have  utilized  the  concept  of  the 
sound  path  in  computing  the  approximate  temperature  in  the  upper  portion  of 
the  stratosphere.  Esclangon*,  in  a publication  for  the  French  Army,  has  studied 
the  problem  of  sound  propagation  along  theoretical  lines.  The  present  inves- 
tigation, though  carried  on  independently  of  Esclangon 's,  parallels  his  work  in 
certain  fundamental  derivations.  However,  none  of  these  investigations  carries 
to  a complete  solution  the  net  effect  of  the  meteorological  conditions  along 
the  entire  sound  path  in  a manner  suitable  for  application  in  sound  ranging. 
This  is  due  in  part  to  the  different  nature  of  the  problem  with  which  they  are 
concerned,  and  in  part  to  assumptions  which  are  not  Justifiable  from  the  stand- 
point of  the  accuracies  desired  in  the  present  sound  ranging  program  in  the 
U.  3.  Army.  We  paramount  concern  0/ tAe  present  invutigation  ic  t/it  evolution 
of  tMoory  of  sound  propagation,  tparing  no  gcnorality  wAsrsver  tAts  endeavor 
i*  compatible  mitb  an  explicit  eolation,  yet  wttAin  yul,  eight  of  the  fact  that 
ite  eventual  utility  can  only  be  that  of  a practical  field  sound  ranging  tech- 
nique. The  solutions  obtained  permit  any  manner  of  variation  with  elevation  of 
temperature,  humidity,  wind  direction,  and  wind  velocity,  these  four  elements 
considered  in  any  combination  whatsoever. 

It  is  shown  in  subsequent  sections  of  this  report  that,  in  long  range 
propagation  of  sound,  the  sound  path  is  somewhat  similar  to  the  trajectory  of 
a projectile.  That  is,  sound  propagation  through  the  atmosphere  is  a three- 
dimensional  problem.  Hence  a study  is  required  of  the  microphone  array  which 
affords  means  of  determining  the  direction  of  sound  approach.  To  this  end 
Section  II  is  devoted  to  a mathematical  analysis  of  the  geometrical  arrangement 
of  microphones  sufficient  to  determine .uniquely  the  space  orientation  of  the 
direction  of  arrival  of  sound. 

*t.A.  tatae.  amed  few*  in  the  Almoephere,  mi.  mo-  *t.  PO.  tv-tta.  July  tssi. 

* 0.  antentert.  Die  SchallauebretUmg  tn  dcr  ttmegHin,  BMObaeb  Mr  OMpSyttk.  Bud  9.  pp.  S9-146. 
s g.  sorIusm.  L'/kouettqee  dee  Cdmone  et  dee  Pnjeetilee,  Henotial  de  VMrtillene  F*mcaiee, 
nr,  89  rue. 
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lh«  solution  of  the  problsa  rsqulrss  slso  s knowlsdgs  of  ths  vsloolty 
i ' of  sound  snd  its  vsrlstlon  with  sisvstlon.  Ssetion  III  includss  s discussion 

I of  squstlons  rslsting  ths  vslooitjr  of  sound  ss  s function  of  tsapsrsturs  snd 

I hunidltf.  Ibsss  nstsorologlosi  slsasnts  can  be  dstsnslnsd  bjr  Mtsorogrsph  ss- 

esnts  In  ths  staosphsrs.  Ths  final  aquations  are  put  In  such  fora  as  to  utilise 
quantities  used  In  current  neteorologlcal  prsctioe. 

The  nonhoeogenelty  of  the  e^lua of  sound  trsnsaisslon,  vis.,  theataos- 
phere.  Is  considered  In  Section  IV.  The  general  equations  of  notion  are  derived 
for  the  propagation  of  sound  In  such  a nedlua,  after  the  nanner  of  developnent 
by  Milne. 

These  equations  of  notion  are  rendered  to  an  explicit  solution  in 
Section  V for  the  case  of  a stratified  atnosphere.  A technique  la  Introduced 
In  the  integration  of  the  aquations  of  notion  which  pemlts  of  parsaetric 
I solutions  for  the  tine  of  travel  of  the  sound,  and  for  the  coordinates,  x,  y, 

s,  of  a point  on  the  sound  wave  front  at  any  position  along  Its  path  of  pro- 
pagation. As  a consequence  the  nethod  of  analysis,  with  but  slight  nodlflca- 
tlon,  lends  Itself  to  the  enploynent  of  graphical  or  nechanlcal  devices,  ob- 
viously advantageous  in  that  laborious  conputstlons  ore  terrlnated  once  tne 
parsaetric  solutions  are  adequately  represented. 

Certain  slnple  cases  of  atnospherlc  stratification  are  considered  In 
Section  VI. 

In  Section  VII  a aathaaatlcal  solution  Is  derived  for  the  position  of 
the  sound  source.  The  solution  leads  to  the  deteralnatlon  of  aparaaeter  which 
aay  be  utilised  as  an  index  of  reliability  for  the  ranging  operations  perforaed. 

No  atteapt  has  been  asde  In  this  prellalnary  report  to  Incorporate 
correction  terns  for  the  tine  of  travel  sad  the  range  and  aslnuth  of  the  source 
of  sound  due.' to  abnomally  high  sound  velocities  In  the  vicinity  of  the  eound 
source,  when  that  source  Is  an  explosion  of  considerable  nagnltude.  Treat- 
aent  of  thle  phenoneoon  Is  deferred  until  such  tine  as  the  essentials  of  the 
theory  herein  are  put  to  rigorous  test. 

Soae  extensions  and  s(>plication8  of  the  basic  theoretical  deductions 
contained  herein,  particularly  in  regard  to  the  tactical  utility  of  neteoro- 
logy  and  forecasting  for  souixl  ranging  and  artillery  fire,  are  tentatively 
pr<^»sed  to  fora  the  sid)staoce  of  Voluae  III  of  this  series  of  reports. 


siie*  m*  for  MMMsias  stMoroloslval  MmKUas  for  mbs  srostsKlos  tsns^ 
Svnsai,  mt  wtu  M SiseasMS  la  9um$  II  of  fUs  aorloo  of  lesorta. 
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M.  TKMT  FH  SPACE  SCmillMTlM  CF  TIE 
• ICECTIM  IF  SUM  PtIPAIATIM. 

I.  Dlrnetlon  CmIiim  of  RotmI  to  Soul  V«v«  Front.'— Consider  a left- 
handed  cocwdinate  ayatea  as  in  Figure  1,  and  an  eleaent  ASC  of  any  plane, 

ehoae  nonal  n has  the  direction  angles 
a,  p.  Y with  the  positive  /-.  T-.  and 
p to  « a i Y Z-Mes  respectively.  If  ee  set  X “ ««*  a 


u • c««  9,  V ■ Y,  then  A,  u,  v becoae 
the  direction  cosines  of  the  noraal  to 
the  plane  iJC.  and  the  equation  of  the 
plane  is  given  by 

X*  ♦ t»y  ♦ V*  = P,  ( 

■tere  p is  tte  distance  along  the  noraal, 
•massing  through  the  origin  P^,  to  the 
plane.  If  non  P^.  P,.  P.  are  any  tfjree 
points  in  space,  and  if  the  direction  of 
propagation  of  the  plane  is  invariant 
ovir  the  diaensions  of  the  tetrahedron 
foraed  by  P,.  P^.  P,.  P,.  tr*n  the  dis- 
tances froa  the  origin,  P^,  p,.  P,  of  ttie 
planes  parallel  to  iJC  through  points 
' 1 , P..  P,.  P,.  are  given  by  equation  (I) 

/ I Fluuts  1.  »jtt,  ttie  coordinates  of  the  given  points 

substituted  for  s,  y,  s.  Boaever,  if  F is  the  velocity  of  propagation  of  the 
plane  ilC  along  the  noraal,  and  if  m let  the  arrival  tiae  of  the  plane  at 
bo  T-.  and  T^,  f,,  F.  at  P-,  P,.  P,  respectively,  then  P.  = f(r^-r^). 
p,=  '(Fi-Fo)*  F,  = *(fa-fe3.  and  ae  have  the  three  equations 


Mwv 


*»x  + y,u  + *»v  = F(r^-r,) 

*,X  + F,U  «,V  = F(r,-F,) 
*,X  + y,u  * S,v  = F(f,-F,) 


We  can  noa  consider  X,  u,  v variable,  though  subject  to  tta  trigonoaetrical 
relation  X*  u*  v*  ■ 1.  and  solve  the  equations  (2)  for  X,  u,  v in  terns  of 


the  coordinates  of  P^.  P,,  P,.  arrival  tines  f,.  T^.  f,.  f,.  and  F. 


The  condition  that  there  exist  a unique  solution  for  X,  u.  v of  equa- 
tions (2)  above  is  that 


If  s.  = s,  = s.  = 0,  then  9 = 0,  so  that  the  condition  that  equations  (2)  have 
a unique  solution  reduces  to  the  restriction  that  any  one  of  the  points  P.. 

Ft>  Fg,  P„  Bust  not  belong  to  the  plane  of  the  other  three.  ° 


U 


SdEfcBfcCuX. 

«lli«  f 0 t ARTILLE 


Awided  timt  are  bo  cbOBon  with  respect 

cooditioa  0 / 0 is  satisfiea,  we  can  write 

Vfo  • 


to  that  the 


and 


A = I 
P 


D 


v=  £ 
P 


V^o 

V^o 


V^o 

V^o 


y. 


J't-rp 

V^O 

V^o 


(*) 


Thus  if  we  can  detenine  f,  we  have  a unioue  solution  for  X,  u.  v. 
can  readily  do  by  aecuis  of  tbs  relation  X*  -f  u*  4-  v*  ■ 1 , or 


This  we 


P" 


V**©  yi  *1 

t 

*1  V**0  *1 

f 

*1  y»  V*"© 

ya  »t 

+ £^ 

ft* 

*t  yp  V*‘o 

^-^0  y. 

u 

'p 

u 

'p  yp  J'p-^'p 

= 1. 


(S) 


Hence  Xi  ii>  v can  be  deterained  uniquely  in  terns  of  determinants  of 
order  three  forwd  by  elenents  contained  in  any  square  array  of  the  system 

*1  *1  V^) 

*t  yt  *•  •• 

*•  ys  V*"© 

so  that  the  direction  of  sound  approach  is  independent  of  the  velocity  f of 
the  wave  front. 


1.  Azlastk  and  Inclination  of  Sasnd  Vave  Front.— It  is  possible  to  convert 
the  direction  cosines  x.  u.  v defining  the  position  in  space  of  the  nomal  to 
the  wave  front,  into  an  elevation  angle  and  an  asiauth  ai^le.  To  this  end  let 
g be  the  angle  between  the  positive  F-axis  and  the  projection  of  the  nomal  PK 

1£ 


on  tte  Xr-plue.  Contldeic  e poaitive  «fasn  aeasurad  froa  tbs  r-axit  toMrda 
tba  poaitive  X-axla.  Tten  froa  Figure  2 ae  see  that 


tM  e = txfM  = 

’ V,  V, 


/ ^ ££1-8  = A . 

*~YvAe  FP«»P  u 

/ ' Hence  6 * tea  x/u,  sod  wo  have  thus  de- 

p g terain^  the  asiauth  of  sound  approach  as 

/ * \ projection  of  tno  nornal  to  the  wave 

l/  \ front  on  the  plane  tangent  to  the  earth  at 

^ * w-X  origin  P^. 

® '^*FICDR1  2. 

Since  v is  the  cosine  of  the  angle  between  ttie  positive  Z-axis  and 
the  poaitive  direction  of  the  nornal  to  the  wave  front,  it  follows  that 
eoi'Vi-v*  is  the  angle  of  inclination  of  sound  approach  with  reference  to  the 
plane  tangent  to  the  earth  at  the  origin  Pg. 

S.  Apflleatioa  to  Field  Operat Ions.— The  above  theoretical  deductions 
“ are  ianediately  applicable  to  sound  ranging  technique  provided  tie  restric- 
tions for  a unique  laathenatical  solution  of  ttie  equations  are  observed-  In 
equations  (2)  above,  we  have  the  four  variables  X,  u,  v,  and  F,  but  by  use  of 
the  trigononetrical  relation  X*  + u*  + v*  ■ J.  we  can  reduce  the  nunber  of 
variables  to  three.  Bence  three  detector  aicrophones  are  required,  one  at  each 
of  the  points  P.,  P.,  F,.  and  an  additional  aicropbone  at  to  determine  the 
time  To  that  the  sound  wave  plane  passed  the  origin.  As  pointed  out  in  sub- 
section (I),  the  positions  of  these  aicrophones  must  be  such  that  D i Ot  in 
other  words,  all  four  aicrophones  aubt  not  lie  in  the  saae  plane. 

In  addition  ttere  are  a few  other  liaiUtions  on  the  placeaent  of  the 
aicrophones.  In  the  analytical  treataent  above  it  was  assuaed  that  the  di- 
rection of  the  aorasil  to  the  plane  reaained  fixed,  and  ttiat  the  velocity  F 
along  the  nornal  reaained  constant,  over  the  diaensions  of  tto  tetrahedron 
foraad  by  Pg,  P.,  P..  F,.  Thus,  in  order  Wjat  the  assuaption  of  a plane  wave 
front  be  Jivtified  and  that  the  velocity  not  vary  over  the  detector  spread, 
the  diasnsions  of  the  space  array  should  be  saall  coapared  with  the  range  over 
which  the  «w>nw<  is  propagated.  On  the  other  hand,  the  detector  aicrophones 
should  be  placed  far  enoi^h  apart  so  that  ttie  tiae  differences  so  secured  will 
be  large  c<»ared  with  the  errors  in  determining  the  actual  tlae  of  arrival  of 
the  wsTO  front.  Also  it  is  desirable,  though  admittedly  not  always  feasible. 


be  as  nearly  ocwarable  as  possible  with  the  horizontal  distances  between  de- 
tectors, for  the  ssae  reasons  as  cited  above.  It  thus  follows  that  for  a 
given  range  and  for  a given  set  of  recording  equipaent  there  is  an  optimm 
gpmc4  mrrmy  whose  diaensions  are  set  by  the  accuracy  of  the  equipment 

in  recording  the  tiae  of  arrivals  and  by  the  accuracy  to  which  the  record  can 
be  read,  and  whose  aaxiaun  dimensions  are  set  by  the  range  over  which  operations 
are  made,  so  that  the  assuaption  of  a plane  wave  and  invariance  of  its  direc- 
tion «d  velocity  over  the  dimensions  of  the  space  array  be  a justifiable  one. 
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III.  Til  TIINITieU  VIIKITT  IP  INII 


«.  Tlw  liloclty  «f  Sound  In  Dry  Air.— It  dm  boon  ohowi  by  thoowtle^ 
ooBaldorntlaas,  nad  oaboUatUtod  ^ oxporiMat,  tbot  tho  voloelty  of  nooBd 
In  n tn  by 


•-YT 


Y = ^ 


whom  n It  tho  prottnro,  p tho  dontlty  of  tho  goo,  tnd  c.  tad  e,  ^ ****  _ 
tpoelfie  httU  of  tho  gM  ot  conoteat  protturo  tad  »t  eontut  voIum.  roip^ 
tivoly.  If  noo  tho  got  it  to  idotl  gtt.  than  Itt  oqottion  of  tttte  It  glron  by 
p m piT , whoro  r it  tht  tbMlato  tonportturo  tad  I it  tho  gtt  contttnt  tppro- 
pritto  to  tht  gtt.  Wt  that  bttrt 


✓V  if  . 


(T) 


Air  it  otooatitlly  t aistiiro  of  gtott,  of  lAich  oxygon  tad  nitron 
aooooat  for  tboat  00f,  tad  tptrt  froa  «tritblo  ottor  vtpor  eoatoat  tad  lucti^tod 
pollatioa  in  iadaotritl  rogioat,  lo  eta  rogtrd  tir  (tt  lottt  in  *bt  low 
phtro)  tt  t naifom  aixtaro  tad  ta  idotl  gtt,  httriag  tho  toM  oqwtiM 
tot  olnglo  idotl  gtt,  tet  t gtt  coaotnat  idiloh  it  t preporly  noightM  tritbaotio 
atta  of  tho  iadividotl  gat  eoeataau  of  tho  gtoOt  praatnt.  Thm  ''•^btot 
Y >^1.408  tad  t • S.870B  X 10*  c.g.t.  uaito  for  dry  air,*  to  that  tht 
of  tt  0*  C.  ta  givon  by  aquation  (r)  abort  it  8S1.87  notoro  par  toeood,  or 

86B.61  yarda  par  ooeoad. 

. Effoct  of  Wttor  Ytpor  on  tho  Voloelty  of  Satnd.— Ut  « bo  tho  partial 
art  of  tho  oator  rtpor;  than  tho  oqaation  of  Itato  for  oator  rtpor  eta  bt 


writtoa 


a 


(I) 


ohoro  a-  ■ Jf.  tho  aolocalar  iiol»t  of 

vppor,  t^  r tho  tboolnto  toaporatoro.  ilUailtrly*  im 

abaolatoly  dry  air  otn  bo  oritton 


• tht  dtnaity  of  ottor 
itloa  of  otaU  of 


uL. 


Pt  *’• 


m 


thoro  .«  « 88.fli,  tfci  oabaerlpt  . roforrlag  to  air  on»lo««l7 


tht 


kwmiditj,  tad  dofiatd  by  tho  raUtiea 


0 • 


h 

♦ 


on 


U 


fl*'-,i4^ '.  i'll  'Ali-'J ‘ J - J • 


»c-e  o^T  — 


IfefMib  U* 


an*  M AvalloiAlan  to  taato  iMfiac 


Sabatltatlng  for  p,  and  from  oquotlons  (•)  and  (•)  m have 


.£E—  a 


• a. 


P,  + Pa  «'  r ^ »d  "a 

B'  r B'  r 


• a. 


• ■■  + Pa  "a  • "w  + (P~*)  "a 


p beia|  the  actual  ataoapheric  pressure  and  equal  to  + a.  We  can  write 
further 


= _a 


■«  * "a 

« p + e -2 4 


and  since  ^ s .8221  , or  ^ approxiaately, 

* 

o = X -y 


8 p ^ 2 g 


(M) 


But 


p = a+Pa  = ^p_r+f^pr 


a.  » 


= ^ S P.  Pa^  • 


and  since  p ■ p«  -f  pw  • 


p a RLdl  [j  + (li  - J) 

"a  ■«  P 

= £1aI  [J  + 0.6q]. 


(12) 


Ibus  for  air  at  teiqierature  T having  a quantity  q of  water  vapor  present  we  can 
UM  tte  saw  equation  of  state  as  for  perfectly  dry  tir  at  a temperature  T', 
where  T>  « T(1  +0.8q).  The  quantity  T'  as  here  defined  is  the  virtual  tca- 
Perstare. 


for  an  ateosphere  not  completely  devoid  of  moisture  we  can  therefore 
write  for  the  velocity  of  sound, 

a = /YJf  a *d.8q)  . 

where  \ is  tbs  proper  ratio  of  the  specific  heats  for  the  aoist  air  mass 
lidersd.  If  the  water  vapor  is  uniforely  distributed  throughout  the  air 


(It) 


IS 


I 


I 


E 


■ui.  Me  can  m-itc 


P 0 


^ 4 ^ 

•a  •• 

where  the  subscripts  i and  • refer  to  dry  atr  and  ivatar  vaaer.  respectively,  at 
tesiperature  T.  The  ratios  Pa/e^,  P,/a,axpress  the  nuaber  of  noles  of  dry 
and  water  vapor,  respectively,  contained  in  unit  voluae  of  ataospbere.  Now 
by  substitution 

• = A 


As  “a  + \ i 

■«  "a 

"d 

A ^ • 


In  the  absence  of  the  value  of  \ at  atacspberic  teaperatures,  wc  can 
use  its  value  1.324  at  locfi  c.  for  the  purpose  of  deteraining  the  coefflelent 
of  q in  (IS) . Thus 


J _ i = a _ J = 0.  30PP 


1. 403  5 


Hence 


a =-.  Aj  (J  4 0.5099q)  . 


Equation  (IS)  expresses  the  velocity  of  sound  as  a function  of  the  constants 
and  R'/m.  for  dry  ur,  and  the  variables  atoAspheric  teaperature  T and  specific 

U J J A.  _ ea*a • 


buaidity  q.  Tie  can  re-write  (IS)  as 


a * Aj  ((J  + 0.6q)  - O.OPOJq] 

V 

= _ Aa  ^(r'  - o.opojqr)  . 


where  V is  the  virtual  teaperature  as  defined  by  (II),  and  -O.OROlqt  toe  cor- 
rection tera  due  to  the  slight  variation  of  Y with  varying  gfuantity  of  water 
vapor  in  the  aoist  air  aass. 

The  specific  humidity  q in  the  ataoq>here  will  e<  asst  be  .04  (this  in  a 
saturated  ataosphere  at  teaperatures  around  10^  F. ),  generally  about  .vif  at 
noraal  suaaer  teaperatures  and  buaidity,  and  a'daeidedly  aaller  value  tba  loner 


J 


the  tcap«r«ture,  prctture  or  huBidity.  The  correction  tera  0.090lqT  on  f'  will 
therefore  be  very  anil.  For  q ^ .04,  the  aost  extreae  cue,  the  correction  it 
i.i*  C.  for  u air  au*  coapletely  taturated  at  J8*  C,  at  tea  level  pretaure. 

Thia  correaponda  to  an  error  of  i^ut  two  feet  per  aecond  in  the  velocity  of 
aound,  or  an  error  of  .18  of  1%.  At  noraal  auaaer  teaperaturea  and  huaiditiea 
the  error  ia  about  0.4*  C. , or  about  0.068  of  1%  in  the  velocity  of  aound.  Thia 
error  ateadily  decreaaea  with  decreaae  of  teaperature,  buaidity  and/or  preaaure. 

It  ia  therefore  obvioua  that,  to  a good  approxiaation,  wecanuaefor  the 
equation  expreaaing  the  velocity  of  sound  as  a function  of  teaperature  and  huaidity, 

a =yVyr  (J  + O.Oq). 

where  Y and  8 have  the  values  given  in  sub-section  (*)  above  ana  T(1  + 0.6q)  is 
the  virtual  teaperature,  a quantity  much  used  in  meteorological  calculations. 

Since  q = r,  where  is  the  saturation  specific  humidity  at  teaiperature  T 
(uid  existing  pressure),  and  r the  relative  huaidity,  «e  can  write 


a = 21.946  yards  par  steond.  (17) 


Equation  (is)  is  the  aoat  practical  expression  to  use  in  connection  with  data 
secured  by  ataosptorlc  soundings,  inasvuch  as  the  auantitios  T,  (a  func- 

tion of  T p only),  and  r are  readily  coaputed. * At  best  the  correction  on 
the  velocity  of  aound  for  the  presence  of  water  vapor  in  the  ataosphere  is  a 
negligible  one.  It  is  to  be  noted,  however,  that  since  the  gradient  of  hiMdity 
in  the  ataosphere  aay  be  of  considerable  nagnitude  in  certain  instances,  the 
additive  effect  of  the  huaidity  gradient  to  the  temperature  gradient  is  to  be 
considered  In  connection  with  the  refractive  properties  of  the  atmosphere. 
Equation  (is)  indicates  that  the  application  of  the  humidity  correction  is  a 
very  simple  procedure.  As  a matter  of  fact,  the  proper  evaluation  of  elevations 
froa  staospheric  soundings  requires  the  virtual  teaperature  correction,  and 
therefore  the  data  required  for  equation  (is)  is  iMadiately  at  ba^.  AtmMpher- 
ic  diarts  are  available  with  the  virtual  temperature  correction  indicated  there- 
cn. 

Fbr  oouvenlence  of  reference,  a chart  of  the  virtual  temperature 
correction,  P.drs....  for  saturated  air  at  given  values  of  temperature  and 
pressure.  Is  given  in  Appendix  A,  page  67  , of  this  report. 
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IV.  EVVATIOIt  9f  MTIOI  F«R  TMC  PIOFMATIOH  OF  lOHaO 

0.  Diff«r«ntiol  EoMtlofiv  of  • Sevnd  Roy.^It  bao  been  sbown  in  tbe  above 
that  tbs  velocity  of  sound  in  still  air  is  essentially  a function  of  teapera- 
ture.  Since  teaperature  aay  vary  over  all  space,  then  in  general  at  a point 
y,  M tbe  velocity  of  sound  is  given  by  a(x,  y,  «).  If  tbe  air  is  in 
■otion  tbere  are  to  be  coiwidered  in  addition  tbe  components  of  toe  velocity  of 
tbe  mediua,  given  by  u(x,  y,  s),  v(x,  y,  s),  w(x,  y,  x).  If  X,  p,  v are  tbe 
direction  cosines  of  the  normal  to  tbe  save  surface  at  Ptx,  y,  s,  then  we  have 
for  tbe  equations  of  aotion,  as  indicated  in  Figure  3, 


^ = X«  + «; 


^ = pa  + w; 


^ = va  + IK 
it 


Equations  (K)  represent,  in  general,  tbe  differential  equations  of  a family 
of  space  curves*  in  tbe  ataospbere.  Ve  can  define  a sound  ray  as  one  member  of 
this  faaily  of  space  curves,  where  each  member  bas  an  appropriate  "ground  value" 
X,,  Pg.  Vp,  and  a unique  variation  of  X,  p,  v along  tbe  ray. 

7,  VariatioA  of  Direction  of  Propaqation  along  the  Sound  Ray.— Tbe  space 
variations  of  a,  u,  v,  and  *>  can  be  determined  by  sufficiently  dense  meteoro- 
logical soundings  of  tbe  ataospbere.  However,  to  integrate  equations  J 1 8)  we 
also  need  to  knM  tbe  variations  of  X,  p,  v with  x,  y,  • ed.ong  tto  ray.  To  this 


/ 


end  consider  two  neighboring  sound  rays,  as  JI3  and  A>B>  in  Figure  4,  an  in- 
finitesimal distance  apart.  Let  55'  represent  the  surface  of  tbe  wave  front, 
and  Ptx,  y,  u andF'is  6a,  y + by • • -F  6*  two  neighboring  points  on  the  wave 
surface  at  time  t.  Then  we  have  for  an  infinitesimal  variation  over  the 
surface, 


*WMW  aecwrauiy  •ums.  Ocm 
eaattasdw  y«*r  tSt  iwflM  tiw 
cat  isaet  aartt  will  as  ■ 
esM.  «t  eallad  • MSMi  tpaat 


tp  tsea 
peiat!  ( 


asaet  eumt  aa  easdlttaa  taat  a tel  s,  w 
•ea  evnra.  MStravar  a aM/ar  14  «l  w art 
It!  aaa  taarafora  aay.  tp  laalwy  mu  taa 


14  t»  art  aaaryaaara 
■a  siaeeatUBMaa. 
w twe-«iaaaaloaal 


stlaa.  aaa  Stawart  1 
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A IkMrr  h^wptti  TkrMfk  tta 


Mrf  an  Aaalicatian  to  SoanA  Kanfini 


<Ci-±^aK}  . X.  + . + 5(x«  + •)  = {f  + 

= uo  + V + 6(u«  + »)  = ^ + ' 

<l(i  t ftil  s VO  -t  w -t-  6(vo  + w)  = ^ + ^(6*)  J 

4t  at  at 

Subtracting  eorrenponding  pairs  of  equations  in  (••)  and  (i#) » we  have 

j;(6»)  = 6(X*  + *•)  = 6*  5j(X«  + •)  + 6>  L(X«  + «)  + 6*  J;(X«  + «) 

^(5y)  = 6(u«  + V)  = bjr  |^(u»  + v)  + 6*  |j(u«  + v)  + 6*  ^(u«  + v) 

^(5s)  s 6(vo  -t  w)  s 6«  ^(vo  +»»)+&«  |-(vo  + *»)  + 6y  ^(vo  + *») 

at  o9  0* 


Differentiating  fbe  relation 


X6»  + u6y  + v6*  = 0 


with  respect  to  tiae, 

j6*+Jtt6y+g6*+X  j^(5x)  + u j^(6y)  + V ij(6.)  = 0. 

Substitutii«  froa  (to)  into  (22) . ee  have 

f^6s-tX6«  f*(Xe  + •)  + X 5y  ^(X«  + «)  + X 6-  + *) 

at  at  • at  9m 

+ u by  + •)  + 4 6»  |j(4«  + v)  + U 6* 

+ v 6»  *-(vo  + e)  + v 6*  f-(vo  + w)  + v 6y  f-(vo  + a.)  = P. 
go'  8«  oy 

is  can  write  (22)  as 

6*  ♦ X 6»  *j(Xs  + ■)  + U 6*  |j(uw  + w)  + V 5*  |j(vo  + w) 

+ ji  by  + X by  ^(xs  ♦ S)  + U by  Jj(uw  + v)  + v by  J^(vo  + w) 

4.  ^ (a  t>  X b«  S-(X*+w)  + U bs  ?;(»*•+*')  ♦ v b»  fr(y*+*'>  = 
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!fe  can  write  2_  (X“  + “)  = a_(X«)  + Sn  , and  similarly  for  corresponding  ex- 
ax  ~ 8x  dx 

pressions  in  (2«)  , so  ttiat 

^6x+6x[\§H+n^-l-v^  + \ §-(\a)  + u |-(wo)  + V (vo)l 
dt  ox  ox  ox  ox  ox  ox 

+ ^6y+6y  [xf*^  + u|2!^+v2K  + x f-(xo)  + U f-(uo)  + v f-(va)l 
dt  dy  dy  dy  dy  dy  dy 

+^6x+6*  tx|^+uf^+vf^+X  |-(Xo)  + U |-(U“)  + V f-(vo)1  = 0. (25) 
dt  dz  dz  dz  dz  dz  dz 

It  is  possible  to  write  A 2_(Xo)  + u 2_(u<»)  +v2_(vo)  more  simply. 

3x  dx  dx 

Consider  a = /(x,  y,  z)  and  A = F{<^:  then  _ 2A  2a  . ffence 

dx  da  dx 

X |-(Xo)  + u |-(ua)  + V |-(vo) 

9x  9x  9x 

= A [o2A2«  + a2S]  + ^ [o2u2a  + u22]  + v [a2i^22  + v2“] 

9a  9x  9x  da  dx  dx'  da  dx  dx' 

dx  da  dx  da  dx  da  dx 

= 22  (j+a(A2A  + y2u  + v2i!)]=22 
9x  ^ ^ da  da  da  dx 

Substituting  in  (25)  , 


dt  dy  dy  dy  dy 


+ 5,[^  + 22  + Aait  + u2H+v2S]=0.  (26) 

dt  dz  dz  dz  dx 

Equation  (26)  must  satisfy  the  relation  (21)  , so  ttiat  ttie  terms  in  the  brackets 
in  (26)  are  proportional  to  A,  u,  v,  respectively;  or 

l[^  + 2H  + A2M  + y22  + v2!!!] 

A dt  dx  ^ dx  dx  dx' 

= i [^  + |a  + A + U + V . (27 

u dt  dy  dy  dy  dy  ' 

= i [^  + 2a  + A 2Ji  + u + V as] 

V ‘dt  9«  ^ 9«  9«  dz' 
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Equations  (27)  give  tt«  variation  of  u,  v along  tt*  sound  ray  in  terms  of 
trie  gradients  of  a and  u,  v,  w,  and  tijeroforo  implicitly  in  terms  of  x,  y,  a. 
TtiO  SIX  equations  (1$)  and  (27)  are  sufficient  to  determine  tlie  position  in 
space  of  all  sound  rays  emanating  from  a given  source  of  sound.  It  follows, 
tt«n,  ttiat  ttie  source  of  sound  can  be  located  as  tt«  intersection  of  any  two 
particular  sound  rays. 

In  Section  II  a mettiod  was  indicated  wtiereby  tt«  boundary  values 
X > u t V.  could  be  determined  for  ttie  paurticular  sound  wave  front  incident 
over  ?fje  detector  spread.  It  is  ttie  purjxjse  of  tt»e  following  sections  to  sug- 
gest a tectmique  wtereby  equations  (H)  and  (27)  can  be  rendered  towards  an 
explicit  solution,  once  tbe  "ground  vfcluos"  X,,  Uo>  known. 


"S  E C R E-^F 


FOR  RRTILLE 


I V.  TIE  FIOPROATiei  OF  SOOlO  II  A STRATIFIED  ATNOSFIERE 

I ‘ It  wai  ■nnlioned  in  Section  IV  that  tbe  source  of  sound  can  be  located. 

[ as  tbe  intersection  of  any  teo  particular  sound  rays.  Tbe  nature  of  equations 

' (>T)  tiOKever,  aakee  it  difficult,  or  at  least  very  tedious,  to  solve  two  sets 

' of  six  equations  eacb,  as  given  by  (II)  and  (27)  . 

A decidedly  sinpler  approach  to  tbe  problem  is  to  trace  eacb  par- 
ticular sound  ray  through  tbe  atmosphere,  making  use  of  the  lam  of  refraction 
and  tbe  variation  of  a,  u,  v,  and  v with  space  {mtttorologieal  corrections), 
and  thereby  determine  tbe  distance  of  the  sound  source  from  tbe  space  array 
and  the  angular  correction  to  be  applied  to  the  observed  azimuth  of  sound 
approach  at  tbe  space  array.  A similar  procedure  for  another  space  array 
suitably  placed  still  determine  a range  and  corrected  azimuth  for  that 
array  for  tbe  same  sound  source*  and  the  intersection  of  the  two  corrected 
. azimuths  will  locate  tbe  sound  source  within  narrow  limits,  provided  tAa| 
due  allotiauce  is  made  for  the  difference  in  elevation  of  the  sound  source  and 
the  origin  of  the  reference  coordinate  system. 

t,  Eqsatlons  of  Notion  In  tho  Stratified  Atmosphere.— Let  us  assume  ttiat 
I as  a first  approximation  the  atmosphere  is  stratified,  or  in  otiier  words  hori- 

zontally homogeneous,  in  the  lower  portion  of  tiie  atmospiiere  through  which  we 
desire  a detailed  an^ysis  of  sound  propagation.  Ttmtn  w = 0 and  u,  v.  and  a 
are  functions  of  s only,  so  tiiat  under  these  conditions 

aus2l'  = 2“  = 22  = 22  = 2s  = o; 
dx  ds  By  By  Bx  By 

and  accordingly  equations  (<R)  and  (27)  above  simplify  appreciEd>ly  to 


^ = ,\a  + u;  ^ = uo  + w;  ^ = vo; 
at  at  <Xt 


\ at  u at  V at  da  da  ds 


(20) 


Integrv’^ing  I iX  m I in  , ee  iiave 

\ dt  li  dt 


or 


log  X = log  u + A 


— = constant. 
U 


Referring  to  equation  (0)  in  Section  II,  it  is  readily  seen  that  the  condition 
x/\S  “ constant  is  simply  the  condition  that  tAe  antnwtA  of  the  normals  to  the 
mmue  front  surface  he  the  same  at  every  point  of  tA«  ray. 

*he  dournlasd  fr«n  e eontonr  nip  la  Um  slelaUir. 
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9.  Th«  Law  of  Refraction.— It  is  ttms  possible  to  further  simplify  our 
equatioM  of  notion  by  choosing  a coordinate  system  such  tliat  for  a particular 
ray  u ■ 0;  or  in  .other  words,  so  that  the  normals  to  ttje  sound  wave  surface  at 
every  point  of  a particular  sound  ray  be  contained  in  any  of  the  family  of 
planes  parallel  to  ttie  XZ-plano.  This  involves  simply  a clockwise  rotation* 
about  the  2-axis  of  our  previous  coordinate  system,  of  amount  ^ + fl  , wtiere 

2 

(?  ■ are  tan  \/|io  • The  direction  cosines  of  ttie  normal  to  ttie  wave  surface 

with  reference  to  the  new  axes  are/l-v’,  0,  v,  since  v remains  unchanged  during 
the  specified  rotation.  Equations  (28)  can  now  be  written 


dx, 

dt 


= /1-v^  a + u^;  ^ 


is.  = vo; 
dt 


(29) 


dt  y dt  de  de 


wtiere  the  subscript  fl  refers  to  coordinates  of,  and  components  along,  ttie  new 
coordinate  system. 


Consider  now  the  last  equation  in  (29)  , This  can  be  written  as 


d_  _ /i-v*  i)i  + /J-v*  da  ^ icy!  du^  . (30) 

dt  y dt  M dz  y dz 


dut 


SO  ttiat 


d_  = 4 . 


dt 


2 /i-\ 


dt 


r -Mrjil  sL  /ir^  . 

dt  V dt 


Substituting  in  (80)  we  obtain 

i.  ^ 

dt  y*  dt  y dz  V 


or 


f /CT  = V g*] 

dt  dz  dz 


(31) 


Mopt  tlM  eOBVtBtlOB  tut  tlM  POtltlV*  I-Ml«  Of  tbO  COOTOlMtO  .JOtU  OfUf  FOtOtK^lnt 
roHlOB  tflMrta  <*leb  tbt  mom  U bolU  propotttoo.  IhU  eoavratlon  luUoo  that 


*»  alMll 

■tii^aiMpp^M’^itivo' iMpwMMt  dr  um  aiai  dr  v. 
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ProB  (io)  we  heve  ^ = vw. 

<lt 


Hence 


IL  /r^  = i^nZ  . 

d#  a dM  dM 


Clearly 


^ /ir^  f/O: 

•'*  VT^ 


_i_^_  a f /ir^  . 

d*  l-v*  dz 


But  from  «2)  we  see  Uiat 


du^  g d 


de  i-v* 


so  Uiat 


d a _ du- 


Integrating  we  have 


/l-vj  /l-v^ 


where  a.  V|  and  w^  are  functions  of  s. 


Equation  «>)  is  the  general  expression  of  the  law  of  refraction  when  both  the 
temperature  and  the  veloci'  of  the  :^ium  are  variable. 


In  Section  II  theie  was  indicated  a means  of  determining  V,  the  tw- 
*tawtM>o«s  vtloeity  of  sound  at  any  Roint  of  tht  sownd  path,  knowing  only  the 
apaue  coordinates  of  the  detector  microphones  and  the  arrival  times  of  the 
sound  at  each  detector.  If  the  normal  has  an  angle  of  inclination  to  the 


horisontal  arc  cos  /J-v*  , tten  clearly  the  apparent  velocity  of  sound  along 


J 


mt  M 4wlloiti«i  t«  $mmi  tm§tm 


tto  tiorisontal  will  be 


Ttie  velocities  V and  a are  thus  related  by  tbe 


so  that  tbe  law  of  refraction  can  be  written 


For  a particular  sound  wave  front  incident  over  a space  array  as  detailed 
in  Section  II,  tfae  "ground  v^ues"  A-,  u^,  v^,  wd  F,  are  re^ily  determined. 
Choosing  the  sound  ray  which  passes  through  the  .origin  of  our  coordinate  sys- 
tem, we  can  apply  equation  (S7)  in  determining  the  inclination  to  the  hori- 
sontal  of  the  normal  to  the  wave  surface  at  every  point  .on  the  ray,  provided 
that  a and  are  known  functions  of  ».  If  there  exists  a level  at  which 

/J-v^  = i,  then  the  sound  wave  is  totally  refracted  at  this  level. 

10.  general  Solution  for  the  Point  of  Total  Refraction.— Lot  us  consider 

in  detail  the  case  where  there  is  ‘a  level  2L..  such  that^  are  eo«  /I-v*  = 0. 
Then  for  the  particular  ray  wo  are  tracing,  we  have  ine  two  boundary  condi- 
tions 

1 

A = /i-Nfe’  y;  u = 0 .*  v = ot  , = o, 
and  A = 1 ; ^ = 0 ; v = 0 ot  •«  = 2^ . 

Letting  /i-v*  = J in  equation  (*T)  we  obtain 

a'=  -fa-  - S . (»•) 

which  is  the  implicit  solution  for  the  level  of  total  refraction  . . 

He  can  now  proceed  to  integrate  our  equations  of  motion*  (*•)  as  follows*. 

- -f  S» 

da ' V v« 


so  that 


- = <ii 


^ l««*l  at  MUl  rafTMtlan  U tM  irwl  «t  Wiles  uw  tatad  rtjr  Smom*  borlnnui  ••  • resalt  of 
SowSlas  of  tao  mmad  iwr  moots  • sorimtal  plaao. 

•ite  MMtive  ot«  la  tat  itfo-lMBd  wwmtr  of  (SS)  la  utroducM  m preotrlr  loeau  tat  paaltiaa  of  tsa 
■tMoaSwrco  la  mm  ar  taa  rafaraaaa  eaordlaaM  apataa  appraprltta  to  ua  aclauts  of  tavad  approaeb. 


I 


I 
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SiBilarly,  - t 


It  it  to  be  noted  ttiat  ttie  integration  with  respect  to  s is  performed  just  up 
to  tbe  level  as  tbe  integrands  in  ttie  expressions,  for  - and-K^  are  not 
defined  for  Tbe  coordinates  (-T^,  Y^,  represent  tl»e  point  of  total 

refraction  in  tbe  atmosphere  of  tbe  sound  wave  surface  whose  normal  lias  an 
azimuth  e,  and  which  is  incident  to  ttie  plane  s = 0 at  (0,  0,  0)  at  an  angle 
arc  cos  v..  It  is  to  be  noted  ttiat  for  sound  rays  incident  to  ttie  plane 
s = 9 at  (0,  0,  0)  at  angles  ottier  ttian  arc  cos  Vg,  the  point  of  total  refrac- 
tion will  not  be  the  same  as  for  arc  cos  v^;  in  general,  each  sound  ray  vill 
have  its  own  characteristic  point  of  total  refraction.  Further  discussion  of 
total  refraction  of  sound  waves  is  left  to  a more  complete  analysis  of  the 
refraction  equation  in  connection  with  audibility  and  zones  of  silence. 

Equations  (39)  and  (*0)  are  not  immediately  integrable,  since  a,  u^. 
and  v^  are  complicated  functions  of  «,  depending  on  the  current  state  of  tlie 
atmosphere  as  regards  temperature,  humidity,  and  wind.  Moreover,  v is  a fuiw- 
tion  of  these  meteorological  elements  and  hence  also  of  «.  It  is  quite  obvious 
that,  even  if  a,  v^,  and  v were  known  functions  of  c,  the  integration  of 
equations  (39)  and  (10)  would  still  be  a hopeless  task  in  ttie  majority  of  cases. 

II,  Generalizations  Regarding  Integration  by  Successive  Layers.—  If  we 
write  our  differential  equations  for  sound  propagation  in  another  form,  in 
particular,  differentiating  with  respect  to  the  cosine  of  the  angle  of  in- 
clination of  the  normal  to  the  wave  front  (which  varies  with  a and  and 
therefore  depends  implicitly  on  «),  we  obtain  a set  of  equations  which  are 
immediately  integrable  provided  only  slight  and  justifiable  assumptions  are 
made.  We  proceed  as  follows: 

From  (31)  we  have 


ijdlZ  = V /T^  [ii  + , 

at  dz  de  ' 


80  ttiat 


a’t  = i.y.lrJtlL . (91) 

If  we  choose  our  interval  of  integration  so  that  ^ and  du^are  constant  for 

dz  dz  ^ 

this  interval,  then  equation  (91)  is  immediately  integrrible.  Let  A represent 
the  lower  boundary,  B tte  upper  boundary,  of  a layer  of  atmosphere  tiaving  these 
properties.  Then  tte  tine  of  travel  of  the  sound  wave  along  the  unique  path  de- 
termined by  tte  angle  of  inclination  of  the  wave  normal  at  the  lower  boundary  of 
ttie  layer  and  the  values  of  temperature,  humidity,  and  wind  ttiroughout  the  layer, 
is  given  by  the  integration  of  equation  (9l)  above,  between  limits  determined  by 
tbe  values  of  the  variable  of  integration  at  t(ie  lower  and  upper  boundaries  of  the 
given  layer.  Eiefore  integrating,  however,  we  need  to  inspect  the  integrand  for 
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■insularities  within  the  interval  of  integration,  and  test  such  improper 
integrals  for  convergence. 

If,  further,  equations  (It)  and  (10)  are  expressed  in  a fora  siailar  to 
(tl)k  the  integration  for  the  coordinates  of  the  point  of  total  refraction 

can  be  accoaplished  by  first  integrating  individually  for  each  layer  of  constant 
da/do  and  du^/d«,  and  then  sunning  the  increnents  of  these  coordinates  for  the 
various  layers  traversed  by  the  sound  wave. 

Let  us  first  simplify  the  appearance  of  equation  (si)  by  the  sub- 
stitution of  new  variables  and  parameters.  -Set  ^ = /J-v', tS«  cottwt  of  tko 
OHglo  of  inclination  of  tht  vav*  front  normal  to  tkt  korioontal.  Furtheraore, 
since  the  interval  of  integration  is  always  chosen  so  that  both  gn  and 

. , . d*  ds 

are  constant,  we  can  set  a = so  that  o is  likewise  a constant  for  any 

interval  of  integration  chosen  in  conformity  with  the  above.  The  term  ^ is 

dm 

the  gradien^of  sound  velocity  in  the  vertical  in  a still  atmosphere,  and  since 
a = ai.94Br'T'  as  given  in  equation  (17)  , it  is  clear  that 

d£  IQ,  in  di:  ^ (St) 

dz  /f7  dz 

where  T'  is  the  virtual  temperature,  expressed  in  degrees  absolute.  The  term 
dT  /da  therefore  combines  the  effect  of  both  the  vertical  temperature  and 
humidity  gradients  on  the  variation  of  sound  velocity  with  height.  Provided 
we  substitute  for  /T'  in  equation  (S2)  a properly  chosen  mean  value  for  the 
layer,  it  is  evident  that  the  condition  * = constant  will  be  very  closely 
approximated  by  a layer  of  atmosphere  with  constant  gradients  of  virtual 
temperature  and  wind  in  the  vertical.  Thus  from  the  standpoint  of  the  meteoro- 
logist the  choice  of  the  intervals  of  integration  is  greatly  simplified;  as  a 
matter  of  fact,  the  present  method  of  evaluating  atmospheric  soundings  furnishes 
all  data,  excepting  wind  data,  requisite  to  an  immediate  determination  of  a t 
and  hence  also  the  intervals  of  integration  for  purposes  of  the  proposed 
sound  ranging  technique  contained  herein. 

We  can  now  write  equation  (m)  as 

“i*  = 7^  -=7^® • 

^ f)(a  + n) 


(It) 


0 at  any  point  of  the 


or 


t = -L.  I ii- 


Ti(a  + ri) 


The  integrand  is  unbounded  for  ti=+i,  n=-a  ^ If 


n - 
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sound  p&th.  then  in  consideration  of  the  assumption  of  horizontal  homogeneity, 

It  is  obvious  that  the  sound  path  must  continue  to  extend  vertically  indefinitely. 
Clearly,  then,  the  condition  n = 0 not  enter  into  the  solution  for  the 
sound  source.  In  view  of  our  choice  of  the  reference  coordinate  system, 
satisfies  the  condition  0 < r\  < 1 always.  The  condition  n - i occurs  only  at 
a level  of  total  refraction,  at  which  the  integrand  of  the  right-hand  member 
of  equation  (as)  becomes  unbounded;  but  on  application  of  the  Cauchy  "u  -test" 
for  the  convergence  of  an  improper  integral  over  a finite  interval  it  can  be 
shown  that  the  integral  in  (as)  converges  absolutely,  over  any  interval  in- 
cluding n * provided  that  at  no  point  of  that  interval  ti  = -a-  The  condi- 
tion ri  - -a  can  only  exist  when  a is  negative  and  lal<  J.  Applying  the  Cauchy 
"u  -teat"  it  can  be  shown  that  if  ri  = -a  st  a point  over  the  interval  of  in- 
tegration, then  the  integral  in  (as)  does  not  converge  absolutely.  The  signi- 
ficance of  the  condition  ^ ° ~a  will  be  discussed  in  a later  section. 


12.  Parametric  Solution  for  the  Time- Increment  of  Travel .—Restricting  our- 
selves to  intervals  at  no  point  of  which  = ~a  , we  can  integrate  equation 
(as).  By  use  of  formula  #228  in  Peirce's  A Short  Table  of  Integrals,  we  can 
expand  the  integrand  in  (as) , as  follows 


— = 1 r — i 

J /J-n’n  (a  + n)  “I  “I  (c 


(a  + r|)  T)  (a  + r^) 

so  that  by  foTraula  il29  for  the  first  term  of  the  right-hand  member  in 
and  formulas  1196  and  #196  for  the  second  term,  we  obtain 


(ta) 


(aa). 


*•- 


da 

dz 


— ^ log  2 J if  |a|<j 

✓7^  n + a 

- i , ,/  Ia|=l 

a r)  * 

i_± 


✓o'-f 


- arc  sin  if  |al>J. 


f (as 

A 


where  represents  the  length  of  time  required  for  the  sound  wave  front  to 

pass  through  the  layer,  from  the  lower  boundary  A to  the  upper  boundary  B, 
along  that  particular  path  defined  in  differential  form  by  equations  (29). 

Prom  the  nature  of  the  solution  for  the  time  of  travel  it  is  clear  that  the 
length  of  time  is  determined  uniquely  given  only  (a)  the  values  of  the  angle 
of  inclination  of  the  wave  front  normal  at  the  boundaries  of  the  layer,  (6) 
the  values  of  the  vertical  gradients  of  virtual  temperature  and  wind  component 
in  the  direction  6 tnroughout  the  layer,  and  (c)  the  mean  virtual  temperature 
of  the  layer. 

13.  Parametric  Solution  for  the  Range- Increment.— Prom  (29)  we  have 


dx- 

dt 


/J-v’  o + y . 
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With  the  aid  of  equation  (si)  we  can  write^ 


yj-V*  <?  -LiU- 




w V + /J-v* 

as  as 

or  nore  conveniently, 

- ^ = .J 

For  the  layer  il9  we  have 


Wt,  ^ + ,,) 

dg 


(%•) 


r **  r ” 

I ^ I /j.,,  ^ (« + 


n) 


(*7) 


Obviously,  iJefore  we  can o™  integration  indicated  in  (%7),  we  need 
to  express  e and  as  functions  of  rj*  These  relations  are  readily  secure  , 
since,  by  choice,  ^ = C = constant,  and  by  equations  (29), 
iz 


i£  = /ir;7  o 

dt  ' 


Thus 


^ = C /I 
dt 


rn 


a 


ii  = C ^1- 


dt 


= C /l- 


n • 


jdjL. 


dn 


n (*  + n) 


(M) 


For  all  intervals  at  no  point  of  which  = *•  « , equation  («8)  integrates 
readily  to 

loy  o = - log  + log  k, 

n 

where  log  k is  the  constant  of  integration.  Clearly,  (*•)  can  be  written 


^3««  roouwu  6.  pu*  (2B)  . 


(*9) 
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nore  compiclly  as 


which  expresses  a as  a function  of  t)  and  of  tne  constants  a for  the  ^iven 
layer  and  k for  the  particular  path  in  that  layer. 

substituting  in  (50)  from  the  refraction  equation  (37)  we  fins 

(—  - (ti  + n)  = * . 

0 

so  that 

T>  a ^ T> 

0 

Substituting  in  (47)  for  a and  frot  (so)  and  (S2)  w"  haye 


P j a !c 

I • I .^(a  + n) 

4 /V'b  ri_ 

**  _ lo  _1_  I dij 

' ’’.fr 


dT,  , (53) 


h. 

providea  the  interval  of  integration  does  not  include  n “•  -oc.  The  first  in- 
tefiral  in  the  right-hani  me.'nber  is  oroportional  to  the  time  of  travel  along 
the  'ound  pith,  and  readily  integrated,  as  indicated  in  (43)  and  ^ 

order  to  integrate  the  second  term  in  the  right-hand  merahjr  of  (53)  with  the 
use  of  Peirce’s  Tables,  we  expand  the  integrand  as  the  sur^  of  three  terms,  as 

follows: 

v'l-T^^ = I [L  /l=ai  _ L /IcJol  - j , 

n (a  + n)*  “ * ^ c(  Ti  + a (n  + a) 

so  that 

jiL  f/T^-dp-  = =i  [i  I dn  - I dr,  - I ‘‘nl 

n (a  + n)*  ^ «J  aj  n + a J (n  + a) 


ti«  ThrMfii  tkt  It— Ml  m IfplicatlM  to  SoMl  RMfiHf 


Peirce's  Tmklas  reduce  these  integrals  to 


^ I ii(a  + r^)*  Of  a 

“•J  '*•  


- log 


- [- 

n T a 

=i  . 1 [_  log  it  /i-a*- 
^ a n 

d«  . • 


n)  /i-T 


.]  + iO=o! 


By  Beans  of  equations  (**>  and  (*§)  we  can  re-*rrHe  equation  (65)  as 


^ a ri  + a 

ds 

= zk  t - -L  /ir^ 
a ^ n 
ds 


• n “ da  ri  + a 

« d» 


Substituting  the  limits  of  r\  for  the  layer  IB,  *re  have  

d»  ® * 


share  (&X.)*,  = - -T.a  • Ev^^ion  (5S)  is  the  parametric  solution  for  AX.  . 

in  terms  of^fhe  parameters  ®ound  ray,  a for  the 

the  particular  path  in  the  layer.  Ha  shall  refer  to  (dX*)**  and  its 

lent  fliven  by  the  right-hand  member  of  (5$),  as  the  rany«-»scr«e«st.  Purther- 

Mrel  SB  Shan  define  eu  the  ratio  of  sound  velwity  Jri^ient  to  sind  vel^ity 

Sradient.  as  the  gtratifieation  paramtter.  Lastly,  se  shall  refer  to  ds^/d* 

as  the  rtfraetiv*  vind  graditnt,  shich  nomenclature  is  an  obvious  one. 
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Equation  (S7)  c&n  be  written  in  another  forn  by  use  of  the  solutions 
given  in  (so)  and  (51),  fron  which  we  have 

n + « n 


and 


4 = (^  - u^)(j  +D) 

M a.  ^ M 


ik  - , + (!g.  _ « ) D 

'’o  '’o 


= la  > + a . 


(59) 


Substituting  in  (S7)  we  obtain 


- = (“•  - - 7-  • 

a da  n 

dM 


(«0) 


Substituting  as  limits  the  values  of  the  variables  in  (60)  at  the  bottom  A and 
top  B of  the  layer,  we  have 


I-  x.i;  = K.  - - a ^ . 


is  V.  is  r\ 

d«  ® d*  * 


(«l) 


But 


•b=  “a  ^ 

* dt 


+ ^ A«. 


»A  de 


Substituting  in  (81)  we  obtain 

o.  + ^ A* 

t-  -f.J!  = <“*A  ^ dH* 

a.  + ^ A« 

* rim 


iA 

dM 


is  »1a 

dM 


(•2) 


dM 
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or  filially,  if  we  set 


K.  Parametric  Solution  for  the  Or  if  t- Increment.— For  the  corr 
change  in  the  coordinate,  (Af^)  , wnich  we  can  appropriately 

A B 

drt/t-tncrement,  we  have®  from  (29)  and  (3  1) 


3*5  footnote  6.  page 


«S  EC  R^frT 


I 


I 


CM*  of  («•). 


dt 


But  fron  (<il)  m cau  write  (»)  as 


«6) 


dv.  = * 


dn*,  (*  + ’’)■ 

dz 


(67) 


lutegratiug  we  fiud 


a + n 


+ a 


(88) 


where  o is  the  constant  of  integration.  Ecuation  (66)  expresses  the  drift 
velocity  as  a function  of  the  change  of  wind  with  altitude,  the  refractive 
wind  u^,  and  the  parameter  a for  '■he  particular  path  in  the  layer  characterized 
by  stratification  parameter  a and  the  constant  J'g/rij  for  the  sound  ray. 


that 


Equation  (18)  is  the  desired  relation  for  substitution  into  (65),  so 


(69) 
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« TkMry  #»  $mmt  Thrmt^  »• 


Mtf  M U $mmt  iMttat 


HW  expMd  the  iBUtrsuid  in  the  lent  tere  of  (••)  in  order  to  InUgreU,  ne 
followe: 

» 1 I -i-  / dq X I- 

✓7=7  n(«  + n)*  VT=7  *J  <n  + «)  + «> 

(J  f-«7(t,  + «)*/2=7  ’ 

if  leMl; 

= X /Lin X r iA i. 

**  / t/1=7  ^ J (r\  + a)  ✓i-t)'  “ 


= -ly  - 


_i_d=5L_  + X /! 

3 a(n  + o)*  Ja  f (^  + 

(J i — ) I- 1 _i ✓i-.n-  , if  I 

«*  J (a  + n)  ® 

^ r- , , if  lal=J. 


a)/7=7  ' 

if  lal=l. 


(ro) 


^IOI• 


/ 


(t^  + e)  /i-Ti* 


✓i-oP 


to^  2 J t nn  - , if  |a|<i  , 

n + a 


- i ^ • ‘f 


— ^ — arc  *»n  i-f-  Oh  , if  lal>J- 


✓o^ 


n + « 


(71) 


Quite  compactly,  vith  the  aid  of  («5)  , ire  can  irrite 


fi 


-ia = =^t-  ioy  it  . 

(n  + a)  ✓T:^ 


(7J) 
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Substitutiug  iu  (ea)  we  fitid  for  laM  1 


[i=|cd(At)*j  - -^(a  log- 


It  is  to  be  noted  thst  the  above  integrations  bold  for  all  intervals  of 
provided  that  at  no  point  of  the  interval  r)  * -n  . 


A HMry  «r  IwaA  Tlirwuh  th«  MM  mi  Awllntlon  to  SmnM  tenfloi 


Kquatloas  (7t)  are  the  p&rMetrlc  solution  for  the  drift-increaent 
for  sound  props(stion  through  the  layer  AB.  As  in  the  esse  of  the  rsnge- 
increaent,  we  can  write  equations  (7s)  for  the  drift-increeent  in  teras 
of  the  aeteorologipsl  eleasnts  at  the  boundaries  of  the  layer.  Thus  we  have 
froa  equations  (si)  and  (es) 


o = V + 


dv  y 

_ 

dw 

“*o  0 


‘o) 


(7*) 


Otilitlng  the  relations  (89)  and  (7«)  on  substituting  in  (7s)  ve  obtain*. for 
the  case  lel^l. 


“o  + j a da 

® * d. 


JL  /I~nl 


, secA"^n  - -r^  •; — ^ j,  n 
a da  du,  j-a  aa  ' 

d« 


secA  ni 
de 


fit  _i_  A . 

dUT  W da  ’I 

de 


(771 


Substituting  the  values  of  the  variables  at  the  boundaries  A and  B of  tr.e 
layer,  we  have 


dv-  a.  dv-  o, 

— 

B d..  O * 


* du^  * 


(-  = V*,  *B  - -o*  ^ - 


dw^  _j_  fli  _ 1a  tlzH 


dii^  J-«*  ^ 


% ■ 
ds  ds 


(79) 


.TPT.S.10O  IB.  1-f  /l-n^ls  ldantle.1  wlUi  secl)-^  n . 
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Eqttfttiaoa  (7#)  »re  the  p&rsnetric  solution  for  the  drift-increeent 
for  sound  props^stlon  through  the  Isyer  AB.  As  In  the  case  of  the  range- 
increaent,  we  can  write  equations  (7S)  for  the  drlft-incresent  in  teres 
of  the  aeteorologipal  elsaents  at  the  boundaries  of  the  layer.  Thus  we  have 
froa  equations  (si)  and  (68) 

0 


(76) 


Utilising  the  relations  (69)  and  (7#)  on  substituting  In  (7s)  « obtain*. for 
the  case  |al^i> 


»*  1 lit  1 i 

+ r_2.  _ U + -0_  seek  fi  - - — 7,  . T1 

'o  • ^ d* 


secH 


da 

dz 


dg 

1 


<^v»  i ^ . 

du.  1-0*  ^ 

dz 


(77) 


Substituting  the  values  of  the  variables  at  the  boundaries  A and  B of  tr.e 
layer,  we  have 

B * ‘^’^O  “b  . . ^ t 

I-  fal,  = *B  dl^  ^ " d«,  * 

dv.  » f.V.  a.  “bw.  + stcirj^\  _/  Ifi.  _u  i.  tA){t  * stch~'‘‘^k) 

d«.  J-o  n„  •>  T:  a d* 


dv,  j /i-nl  _ ii 

d«^  J-o*  do  ’Ig  do  ’I 

do  dz 


1 . 


(78J 


.WT.M10D  /CTl«  Identical  with  sech'  r,  . 

n 
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* and  bjr  vritlng  the  valuee  of  the  a,  and  v.  at  the  boundary  B In  terns  of 
their  values  at  the  boundary  A and  their  gradient  between  the  two  boundaries, 
we  obtain 


dv. 


r « 1"  < “•’»  . » ^ dv.  **  ■ j_  — ow.  a. 

I-  = (V.,  + &•)  *,  5^*—  t,  + ^ t. 


+ ^ l(^  - "a^  - ^ A.  + + 11£^) 

J-a*  n ^ ® ° da 


(79) 


- ^5^  - “e  + 7-4 

d< 


d£ 

d« 


n,  do  ’1. 

" d. 


But  „ _ da/ds  _ , so  that 


du./d« 


_a._„  +^  = ji_u  + — -S — L. 

♦a  “ »1„  *A 


= -a  - 

n. 


'B 

IS 


'B~“A 


(80) 


Hence 


[-  rj;  = V.,  (*B-*a)  + ^^**8  ^*'’**^  ■ ^ 


t. 


du, 

dv, 


«*  ‘*b“»a~‘*a"»_J)  [(t,-t^)  + -i-(s«cA-^TiB_seeA-^^A)l 

'h  o.-o.  — 

• ■‘-a  'o  “B  A 


1_  [_fL  A,  h.  '-4^1 

d«,  1-a*  «B-«A  “b-“a 


dz 

A 


(81) 
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J 


srtUug 

(AT.)  ^ (i  - 


•■-•a 


At 


_ **1‘>1  f ■ *»~*A  . , (i  - !!!l!4l!£l!i)(,Ac*-‘n»-M«*-‘»u)  (u) 


tM  (eo»"**li)  + — ton  (eof"*nA)]  , 


B-*A 


«*•«*•  l«l  = I *B~*A  I ^ 1. 


•»  *A 


Tollov^ng  the  sane^procedure  for  the^case  lal=lf  *®  fii*i 

. y a /«  + 4l_  a)t  _ 1 2(1  + a)t  - i ^ 4(J  + 9)4-[*«e*"S  + 

- r*  - ^ T 3iiu^v>  ^ (n+«) 


ds 


ds 

dv. 


dM 

dv^ 


(8S) 


-vt-1^2t-lfl!2  2t-li!2a  iifiiLia  _ I filt  2 usA_iu 

• 3 d«.  *1  3 <4*.  ® 3 <iu-  n 4a  3 *1*  * 42 

• • 3*  3a 


. I f!!»  a.  J.  ^£2 
3 3ua  3g  r|  + a 

3« 


Sttbstitating  the  Halts, 

t-  v;  - (...  ^ ^ *•)..  - .....  - 1 ^ (i  - ^ ^)..  *5  ^ (^...1., 


_ 1 £l2  (la-s.  -ii2A*)ii£4IiDl  + 1 (la  _«  )ii£4Il!lA 

3 3«^  '’I,  *A  3*  4fl  ^ '^o  **  ^ 

3fl  A *** 

f!j2  ^ 3s  * + i fi 

-■  A*  _f_  A_I_M  _I_ 


3 3« 


4a 

3> 


3 3*  « 


4a 

3* 


■ 1 5^  i n;~5 nrr; — ’ 

3* 
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dv. 


I-  r,[  • - i - i 

dir 

3 d»^  *A  oj-Oj^  da 


d» 


j dv^  j tan(co5~*^B)  tan(cos~*nA)^ 

3*  du^  ^ M -r  a ti  + o 

d*  ® * 


(8<0 


Hence 

(&r») 


AB 


•a  3 -U^  n„  *A^  3 *B  *A 


■HAt) 


AB 


I *A 


3 -u. 


[^i- 

% 


Bj  ^secA~^^B-5ecA~^^A 


) 


1 [ °B  tgn(cos~*'^B)  _ 

’b"“a  ''b  '*’  ® “B~“A 


tan ( cos 


— ^]1a«  # 


® '^g„-g.  n_  + a “b~“a  + « "J 


(85) 


where  )a|  = 


^ . o.-a 


B “A 


1. 


15.  The  "Lens-Effect"  of  Critical  Atmosphere  Layers.— Let  us  investigate 
the  nature  of  sound  propagation  if  n = - a (where  |®|-J  since  0-r^l  always)  at 
some  point  of  the  interval  of  variation  of  ^ . It  has  already  been  shown  in 
the  analysis  of  the  integrand  in  equation  (H3),  which  represents  the  time  of 
travel  of  the  sound  through  an  atmospheric  layer  along  a particular  path,  that 
the  integral  does  not  converge  absolutely  if  n = - ® at  any  point  of  the  inter- 
val of  integration.  If  n = - « is  an  end-point  of  the  interval  of  integration, 
then  since  the  integrand  is  of  constant  sign,  the  Cauchy  "U  -test"  indicates 
that  the  integral  diverges.  The  significance  of  this  condition  is  demonstrated 
by  the  following  analysis. 


Prom  equation  (32)  we  have 


du„ 


da  = ’J  (dfl  + n _ 
dz  a dz  dz 


) = 2(a  + n) 

dz  ® 


Therefore,  it  is  at  once  obvious  that  -g  is  the  limiting  value  of  , since 
dt^/dz  = 0 when  ^ = - a • 

Suppose,  now,  that  =-a  <J  at  the  end  point  of  the  interval  of 


Integration  n to  n of  equation  (»3) 

A B 
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S’  C C n c T 


j 


(A*)^,  = =^,(*0f  + -ri-  log  2 


• T r.  1 _ ^ «ni  " /i-a*  /J-nJ, 

= :f-  [iof  r ^ iof  2 ■-^^- !^] 


fi^  + a 


■“'{ 


± [log  I*  ./i?  + 
da  ^ 

da 


= ± [log  + 

da  M. 


n + a 


1 •*•  aiT\^~  /i-g* 
n + a 

'i 


_ JL  log  J+  /IrP 

da  “® 


Itm  [log  2 ^-i- 


^1-ar  Tv*-a 


lim  [log  2 

n*-a 


= log  2 + lim  log  '*' 


Differentiating  both  numerator  and  denominator  with  respect  to  ’i  in  the  last 
terra -of  equation  (S7)  in  order  to  evaluate  the  indeterminate  form  0/0  , we 
have 

lira  log  [i— t .aa.=  /i-g 


= Itra  log  (- 


Hence  (At).,  diverges  to  infinity  as  n approaches  -a  . 


But  we  are  reckoning  with  only  finite  intervals  of  travel  time  along 
sound  paths;  hence,  for  an  atmospheric  layer  of  stratification  parameter  a.  • 
we  can  conclude  that  •/  rj  = ■>  ® at  any  roint  of  the  layer,  then  = - a 


S "E'O'^R^E-  T 


I «v«ry  foint  of  that  layer.  On  the  other  hand  if  t>  ^ - o at  anv  point  of  the  layer, 

^ then  n can  never  become  -a  in  finite  time  in  that  particular  layer. 

It  follows  further  that  since  n cemnot  pass  through  the  value  -a, 
and  since  0 < 1 atnospheric  layers  of  stratification  parameters  0 < ~ oi  ^ 1 

I mill  reduce  the  divergence  of  the  bundle  of  sound  rays  about  n “ - a in  the  case 

' da/da  > 0,  and  increase  the  divergence  in  the  case  da/de  < 0.  But  it  is  to 

be  noted  that  this  pbenomenon  is  restricted  to  those  particular  rays  belonging 
to  a family  of  azimuth  direction  0,  where  du^/da  is  of  a magnitude  consistent 
with  0 < - a < 1.  Since  the  refractive  wind  gradient  du^/da  will  necessarily 
vau'y  with  6,  it  is  at  once  apparent  that  for  layers  of  sufficient  vertical 
extent,  the  normal  divergence  of  sound  rays  may  be  disturbed  locally  ed)out  that 
ray  for  which  n = with  a consequent  variation  in  sound  intensity.  The 
discussion  of  the  importance  of  this  analytical  deduction  is  reserved  for  a 
, subsequent  report,  but  it  may  be  pointed  out  in  passing  that  the  above  analy- 

sis proves  a useful  tool  in  the  interpretation  of  the  abnormail  propagation  of 
sound  to  great  distances  and  the  phenomena  of  zones  of  silence  auid  zones  of 
; marked  audibility  at  considerable  distances  from  the  sound  source. 

r I The  solution  for  the  time-, range-  and  drift-increments  for  the  case 

[ I . - a will  be  discussed  in  Section  VI. 

16.  Alternate  Expressions  for  the  Time -Increment. — The  time-increment  of 
f j . travel  for  a sound  wave  through  a layer  of  atmosphere  is  given  by  equations 

[ j . (*5)  in  sub-section  (I2).  The  three  solutions  given  depend  on  the  range  of  a . 

so  that  only  real  variables  appear  in  the  expressions. 

The  time-increment,  nowever.  can  be  expressed  in  another  form  by  change 
of  the  variable  of  integration,  or  for  that  matter,  by  transformation  of  the 
expressions  in  equations  (^5) . Thus 


since  d(cos~‘ri)  = -cIt^  . Utilizing  formula  #839  in  Peirce's  Tables,  and  the 
y'T^ 

substitution 


s = tan  1/2{cos~^t])  , (90) 


k HMry  «f 


iti«i  Tlirnli  tti 


««  M Awlleittw  to  Umd  Iwtlti 


latrodneloi  tte  irrlt* 

<*  - '*'  '*  * S/^” 

Mhlcb  by  use  of  fonatlo  ffl61  In  Polrco'o  Mia*  bMOMo 

fliuA'tiO. , _i_  [i  p^)  - » 1 

/ « + ^ ^ V“  V 

' • 

fro«  trlgonowtry  we  have 


* = *«i  l/2(eo*~\)  » 


80  that  (tt>  becoaes 


= -sit- 


If  |sl=l>  ao  can  aaks  use  of  the  stapler  integrals  AQ6  and  fS07  In  Peirce's 
Tsktst,  securing 


JHtgr*i»i  = JL  di^ 
a + n ® n + « 

..  lol-t 


•'for  |«|=i 

Peirce’s  MU*,  foraula  fSOO,  give  other  variations  of  the  expressions  in 
(n)  above. 

By  writing  ^ p 

/ rfn-  = -/rfCffiTinl  . 

Jr^^  J " 


w«  obtain  by  Peirce's  MU*,  f808. 


= -for  t«.(5  + . 


J 


s F r R^  T 
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H*ne«  thtfr*  *r«  » wld#  varitiy  of  forn#  for  writing  th«  int«gr*tion»  indicatad 
in  oquation  (at)  above.  For  reference  there  are  listed  below  a nueber  of  tljeae 
expresaiona  froe  which  selection  can  be  aade  nost  convenient  for  the  case  ianai»: 


i 1^5?  , if  |al»l. 
* n + « 


H a 


—i~—  log 


n + a 


- log  i-ViriT 


1...  tank-*’  raai- 

i + «n 


(At)^j=  'M  - S«C*"*T 


+ -a2_  /Li-fl 


-log  + ton“*/  ^ >2  ) 

'■#  ^ / J + n " 


-=i—  i -t  ttO 

n + a 


rf-  a 


+ J/ J + 


1 + an 


if  laMl. 


I 

I 


r 


VI.  SOME  SIMPLE  CASES  OF  ATNOSPNEIIC  STIAT IF ICAT 101. 

In  the  previous  section  there  were  derived  expressions  for  the  tlae-, 
rsnge-,  snd  drift-  incre'xents  for  the  propsgstion  of  sound  through  % horitoh- 
tslly  homogeneous  Isyer  of  staosphere,  c^srscterl^ed  by  s strstificstion 
psrsmeter  besring  the  rstio  of  the  sound  velocity  grsdient  to  the  refrsctive 
wind  grsdient  within  the  Isyer.  The  solutions  derived  were  for  the  aost 
genersl  csses,  with  da/dt,  du^/dt,  dv^/di  constsnts  not  necesssrily  zero,  snd 

the  sound  psth  s spsce  curve,  hsving  s point  of  totsl  refrsction  in  the 
straosphere  st 

In  the  present  section  the  solutions  for  the  increments  in  time, 
rsnge,  snd  drift  will  be  derived  for  some  psrticulsrly  simple  straospheric 
strsiif icstions.  They  msy  be  considered  ss  degenerste  c'sses  of  the  genersl 
solutions  of  the  preceding  section;  snd  sre  given  here  inssmuch  ss  the  ex- 
pressions of  the  genersl  solution  sre  inoeterminsie  for  some  of  the  psrticulsr 
csses  cites  here. 

17.  The  Csse  of  n “ ~ * Throughout  the  Atmospheric  Leyer.— Since  ti  is 
constsnt  for  the  Isyer,  tr.e  integrstion  csnnot  be  performed  with  respect  to  ^ 
ss  in  the  more  genersl  csses  shove,  .--owever,  we  csn  perform  the  required 
integration  resdily  by  using  * ss  the  vsnsoie  of  integration. 

From  equations  (29)  we  have 


S T 


«ta*r«  As  * 


A« 


4 7^  ‘"  '■*  * 'J'-**” 

X -i-  i„  2i 

*‘ 


_ 


for 


•■-•a  /I^  -a 
Pbr  the  r«age-lncr«aent,  we  h&ve  froa  equ&tioiw  (tt) 


ctt 


s • + ••  • 


so  thst 


- = J*(n«  + ••)<** 


It  is  required  thst  s snd  be  expressed  ss  functions  of  s in  the  integrsnd 
in  (I0«).  For  we  hsve 

" "♦a  • 

iftifre  is  the  refrsctive  wind  st  the  lower  boundsry  of  the  Isyer.  Froe 
tnis  relstioo  sod  equstion  we  hsve,  on  substitution  in  (lOO) 

yi=^J  a<— a)i 

+ , /* + ^ f- UziAli* 

^ I *i  ^ ^ J ^ £^*-**^1 


MtatM  J.  (SS)  ASM  »t  sew  fW  *✓*  = 0.  «>»•  «“•  »•  “IT*"  “ S'* 

•l»M  q«s^  *s»  *•  tv8ix»aw«  *S*t  U.  it  dc/d>  - 0,  tsw  « • 0,  pro»W«a  dv  / . 


4« 


S'TC"R 


(100) 


(101) 


J 


•Mfhw 


II  -I— — -r- irr-  TTIi  rW 'll'— ■.-■». 


M «M>llQllt«l  to  I 


batwaaa  ih*  Halts  of  tte  lomr  sad  oppsr  booadsrf  obtala 


+ / (*-a)  , 

K ♦g(— a)1 


so  tiMt 


• S*  -3ir;i=;  ’**  * ■ *•  * S'—**' 

«»■  (ii)  /t:;^  *• 


Aa  4^  (At) I 


d»  (da)*  /Ti; 
dz 


l(ag-o*)  - a*  lor  ^1  . 


By  use  of  the  reUtlons  derived  in  (S9)  ve  can  ivrite  (lOS)  aore  coapsctly  ss 


<«.).. ' - 1 


l(At)*8  + ta«(cos'*r))  A*. 


This  result  is  snslsgous  to  the  expression  in  equation  (Ss) . 

The  solution  for  the  drift-increaent  is  siailsr  to  that  for  the  ranje- 
Increaent,  so  that  we  can  iaaedlstely  write 

- («.)„  = * -5^  5?  ^ "••-•*>  - “ il' 

iTi  n 


= (v.  - a.  — t.  ' A)(At)^8  + — I ■ — .. _— 

•a  * a, -a*  "e,-"»A  « 


I 
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Hence 


This  result  is  anslogous  to  the  expression  for  the  drift-increment  in  equation 

(•2)  • 

I8>  The  Case  of  Constant  Rsfraetiva  Wind  and  Variable  Temperature  Throughout 
the  Atmospheric  Layer.— In  this  type  of  stratification,  du^/in  s 0 and  da/da jtO, 
so  that  a-KD  . Hence  for  the  tine  increment  we  have  from  equations  (90), 

- seeh~*^rfj^),  (1C 

“b"“* 

Por  the  range-increment,  we  obtain  from  equation  (6S)  , 


(A/.)..  =-w«(At)..  + ( B-tan(cos"‘T)«) A_ta«(eos"*n.  )lA*. 

(108) 

Por  the  drift-increment  we  have,  from  equations  (29)  and  (31), 


Expressing  as  a function  of  ri  , 


dv^  dv^ 
dt  da 


“i 

• a 


^ - «•)  n 


Thus,  since  dt  = 


n •'■'-n 


and  integrating  we  obtain 


dv^  = - “*)  ‘*’1  , 

0 

dz 


= 5^  a + <y  , 

* ia. 


J 


A ThMrjr  «#  SontiA  rrawt*ti«n  ThroMt'k  tht  AUmpAirt  mA  •»  AMltc«tl«n  to  Sowitf  lo««log 


wfasM  e it  tb«  constant  of  integration,  and  given  bf 


o - 


•■-•a 


Substituting  in  (in)  , 


- - “•>  r * ’ r'- 

J <F'*  J-^  J 


so  that 


f.J*  = — I— A i — A*_  (-00*“*''.  + eoA-^n^]  + A I ^(At)4g 

•'a  0,-0^  »1  S,-0*  " • S,-S^ 


(lu) 


(ll«) 


(Mb) 


Hence 


(&f.)  . = -t  ---a"  “ "o)  l«o.-‘f,|-«**"NAl  f 

AB  ^B*^A  '•i*“*A'  *0 

This  term  is  analogous  to  that  given  in  (|2).  but  not  strictlr  sinilar,  since 
(•2)  contains  the  tera  (o^  ) and  « , which  in  this  cass  approach  sero  and 

infinity,  respectively. 

19.  The  Case  of  Constant  Teaperatwre  and  Vsrisbie  Refractive  Wind  Thros«he«t 
the  Ataospheric  Layer. <»In  this  case,  da/dt  = 0 and  do^/ds  / 0 » Equations  (M) 
cannot  be  used  in  their  fora  to  deteraine  the  tiae-increaent,  since  they  lead  to 
an  indeteminsnt.  But  conveniently,  equation  (»3)  can  be  written 


Hence  for  the  layer  40  , 
(At)is  - : 


_ ^ idS  , 

ds 


[tan(cos‘"^nt)  - tos(coi"*inji)] 


(HT) 


For  the  range-increaent  we  have  froa  (S7) 


(Mi) 


(Mi) 
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20.  Th«  Ca«a  of  Constant  Taoparaturo  and  Constant  RafractJve  Wind  Throo«hout 
the  Ataoaphorie  Layer.— In  this  case  both  d«/d«  and  du^/d*  are  zero, 
a has  no  Meaning.  But  is  constant,  so  that  the  integrations  can  readily  be 
perforMed  with  respect  to  m.  Prom  equations  (29)  we  have 

/■'  ' 

so  that  for  the  layer  AB 

(AO*b  = — ^ 


For  the  range-increment  we  have 


fdx^  = f(r^ 


|0  + 


- ■ Q-~  Az  + u^(At)^g . 

yi-r^ 


(12 


Hence 


(AX^)  = - - cot(cos-^Ti) 

AB 


(13 


For  the  drift- increment, 


dz 


/j-n' 


need  to  express  as  a function  of  z,  and  this  we  can  do  simply  since 


dv^ 


= *'•>  ^ df  • 


Substituting  in  (i3l)  and  integrating, 

I Az 


'•■■I 


Ills 

,4,? 

O /1-ri" 


- I'.l  = I 


= lAsi!  — . 

dz  2 a /T^ 


(I 


o /i-n* 
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*IU  « NATICMTICAL  MMITIOI  rOI  THC  SSttlO  SOWCC. 

Ib  the  foregoing  disciualon  a technique  haa  been  deviaed  vhereby  a 
gltm  xniM  traw  can  be  traced  through  a atratified  ateoaphere  along  a parti> 

^ nanlfold  variation  of  the  Mteorological  elenenxa. 
w conblnation  of  vertical  gradienta  of  tenperature, 

•Ind,  provided  only  that  the  gradient  of  each  neteorological 
eleeent  be  conatant  within  the  atratun  of  ateoaphere  conaidered.  The  aound 
ray  thua  deterelned  ia  a contlnuoua  apace  curve.  Expreaalons  have  been  dsduced 
for  the  tine  of  travel  of  the  aound  along  the  path  traveraed,  and  for  the  In- 
creeenta  in  the  apace  coordlnatea  of  a given  point  on  the  wave  front  aa  a reault 
of  the  propagation.  All  poaaible  inclinationa  of  the  wave  normal  have  been 
conaidered  in  the  analyaia,  up  to  the  level  at  which  the  aound  ia  totally  refracted. 

However,  it  is  conceivable  that  on  occasions  a sound  wave  may  be  in- 
cident over  a locality  without  having  undergone  total  refraction  at  some  point 
along  the  path  of  propagation.  An  obvious  example  is  the  case  where  the  sound 
propagation  is  along  the  horizontal;  another,  where  the  source  of  sound  is 
favorably  elevated  with  respect  to  the  locality  over  which  the  sound  wave  is 
incident.  An  important  consideration  not  to  be  overlooked  la  that  there  may  be 
several  distinct  sound  wave  surfaces  Incident  over  the  same  locality,  but  to 
each  and  every  one  is  attributed  a unique  aound  path.  Using  the  time  of  first 
arrival  at  each  microphone  of  the  space  array  as  representative  of  the  arrival 
of  the  identical  wave  pulse  at  all  microphones,  the  path  determined  from  these 
arrival  times  and  microphone  positions,  and  from  the  existing  meteorological 
conditions,  ia  tht  path  of  least  time  of  travel. 

It  is  the  purpose  of  the  present  section  to  apply  the  rigor  of  mathe- 
matics to  the  problem  of  I'ocating  the  source  of  sound,  given  two  sound  rays 
determined  by  use  of  the  technique  developed  in  the  previous  sections.  It  nay 
be  pointed  out  at  this  Juncture  that  mention  was  made  at  the  outset  in  Section  V 
that  the  'sound  source  can  be  located  as  the  intersection  of  any  two  particular 
sound  rays,  involving  the  solution  of  twelve  simultaneous  equations,  at  beat  a 
tedious  if  not  impossible  task.  Yet  with  recourse  to  elementary  space  geometry 
we  can  effect  a nthematical  solution  for  the  sound  source  in  the  case  of  the 
stratified  atmosphere,  prwlded,  however,  that  for  sound  waves  diffracted  by  a 
physical  body,  the  direction  of  the  normal  to  the  undieturbed  wave  surface  be 
known  before  incidence  on  the  diffracting  body.  The  inclusion  of  this  reservation 
in  connection  with  the  solution  of  the  eight  simultaneous  equatidns  required 
(equations  (2*)  above),  is  an  obvious  necessity,  since  any  phenomenon  which 
diverts  the  direction  of  the  normal  to  the  wave  front  from  that  computed  by 
means  of  the  refraction  equation  will  necessarily  introduce  a cumulative  error 
for  the  remainder  of  the  sound  path  conaidered  to  be  in  the  unobstructed 
atmosphere. 

Ve  shall  therefore  limit  our  solution  to  those  cases  where  the  sound 
path  at  any  instant  is  at  least  800  to  100  fast  distant  (depending  on  the  wave 
length  of  the  sound),  from  any  object  capable  of  producing  diffraction  of  the 
sott^  wave.  Thus,  with  sets  of  equations  as  (29)  , one  set  each  for  two  apace 
arrays  of  microphones,  we  can  determine  the  position  in  space  of  the  source  of 


Moad  vitb  an  nnetaMn  dnpendtnt  on  (•)  th*  accuracy  to  vhlch  the  Inatrunental 
eoulpwnt  will  Indicate  and  record  the  existent  conditions.  (6)  the  degw  to  • 
which  the  ssswiption  of  horlsontal  hoaogenelty  Is  Justlfl^  In  the  portion  of 
Um  ftUKNiDtor#  coQsidsrtfdt  (c)  th®  ®xt®Qt  th®t  th®  th®or®tic®l  tr®itB(int 
herein  presented  represents  the  actual  propagati<»  of  sound. 


ai.  Trsesforestlon  of  Coerdlsetrs.— Consider  two  apace  arrays  of  alcrophones, 
and  choose  one  vertex  of  each  as  the  origin  of  a left-handed  coordinate  syatee, 
with  the  r-axls  directed  northward.  X-axls  eastward.  (See  Figure  W. 
the  rena\nlng  vertices  of  each  tetrahedron  la  space  coordinates  referretl  to  the 
coordinate  systea  selected  for  each  tetrahedron.  For  the  purpose  of  ready 


Identification,  let  the  generalised  coordinates  appropriate  to  space  arry  (1) 
be  *,  y.  *;  those  appropriate  to  array  (2)  be  y».  the 

the  origin  of  the  y»,  *'  coordinate  systea  be  known  with  respect  to  the 
*.  y,  s. systea:  In  particular,  let  P be  the  distance  in  the  ^irection  of  the 
origin  of  the  %•.  y'.  *'  systea,  let  M be  the  corresponding  dlstwce  in  the 
y-dlrectlon.  and  let  E be  the  difference  in  election,  sign  considered. 


Further,  let  e and  e'  be  the  asiauth  of  sound  approach,  deterained  as  outlined 
in  Section  II,  for  the  corresponding  space  arrays. 


The  solutions  for  the  tine-  , range-  , and  drift-increments  derived  in 
Section  V are  expressed  relative  to  a coordinate  systea  . 

orientation  such  that  is  directed  with  an  aziauth  and  sense  the  same  as  tne 
sound  wave  noraal.  The  **  • y*  . **  systea  represents  a clockwise  rotation 

of  aaount  5 + 6 ,of  the  x,  y,  m coordinate  systea.  For  the  other  space  array 
2 

a rotation  of  aaount  5 0'  bas  been  perforaed,  and  the  new  axes  bear  the 

2 

coordinates  x'^,  , y\,  , . 


F 


I 


L 


The  Bound  nve  which  arrives  with  atinuth  e • elevation  eo$  /J-v' 
at  the  space  array  (1)  Ium  a particular  sound  path  appropriate  to  it  only;  like- 
wise for  the  sound  wave  incident  at  space  array  (g).  If  we  consider  indivMual 
strata  of  ataospiiore,  wherein  the  vertical  gradients  of  virtual  teeperature  and 
of  wind  coeponents  for  array  (1),  «.>  . for  array  (8),  are  constant, 

we  can  deduce  the  range-  and  drift-  increnents  tor  each  stratue  by  recourse  to  the 
appropriate,  solutions  for  and  derived  in  Section  V.  To  this  end  we 
also  require  the  tiee-increaent  At  , the  tine  of  travel  of  the  sound  wave  through 
the  layer  along  tht  appropriate  patA.  Inspection  of  the  solutions  for  these  in- 
crenents  reveals  that  they  are  functions  of  various  constants  for  the  layer 
considered,  and  of  the  cosine  of  the  angle  of  inclination  of  the  wave  nomal 
at  t*«  koundarios  of  tkt  laytr  only  , which  we  can  readily  determine  by  neans 
of  the  refraction  equation  (S7)  . And  therefore,  for  a given  layer  whose 
aeteorological  structure  is  known,  and  for  a given  ray  along  which  the  variation 
of  Tf  is  known,  m eon  dttormint  tht  position  in  space  of  the  intersection  of  the 
sound  ray  witA  Aorisontal  planes  representing  the  boundaries  of  the  varioii* 
layers. traversed  Ay  tA*  sound  vave.  The  coordinates  of  these  points  will  be  the 
cuiullative  suns  of  range-  and  drift-  increments,  calculated  from  the  origin  of 
the  reference  coordinate  system,  at  the  appropriate  levels. 

Since  the  points  along  the  sound  ray  appropriate  to  array  (1)  are 
expressed  in  terns  of  , y^  , s^  , and  those  pertaining  to  array  (2)  in 

’ *o>  • required,  that  a common  coordinate  system  be  established 

for  ease  of  reference,  and  for  the  determination  of  the  point  of  intersection 
of  the  two  sound  rays.  For  the  sake  of  convenience  in  position-locating  on 
prepared  naps,  we  shall  use  as  the  standard  system  the  x,y,s  coordinate 
system  already  chosen  for  tetrahedral  array  (1).  The  points  expressed  in  x^  , 
y^  , coordinates  may  be  converted  into  x.  y,  w coordinates  by  simple 
rotation.  Thus 


X = X.  cos  (5+  e)  + y*  *»"  (^+  0) 

2 2 (I3») 

y = y*  CO*  (IJ+  e)  - x^  s»n  (j  + e) 


But  cos  (5  4.  o)  * - »»"  6 , *»«  C*  + 0)  = cos  0,  so  that  equation  (I3H)  can  be 
2 2 
written  more  simply  as 


X = - x^  sin  0 + y^  cos  0 
y = - x^  cos  0 - y^  sin  e 


(136) 


s = s. 


Coordinates  in  the  x^,,  y^,,  s^,  system  when  expressed  in  terms  of  the 
X,  y,  m system  require  both  a rotation  and  translation  for  conversion. 


56 


's-ce^ert 


1 


• S«Mtf  frotiptiw  Tkroaik  tta  M— plan  m<  m t«  %mmi  iMiiaf 


Thus  initially 


But 


Hence 


X = Xt  ^ P:  y = t'  * M:  • = *'+£. 

X'  ~ x>^,  cos  (5  + e*)+  y>^,  stR  («  -f  9') 
2 2 

y'  = y\i.eot  (5  + e')-  x\,  *in  («  + e») 
2 2 

X = P - x'^i  StR  0'  + y'^,  cos  0' 
y - M - x'^i  cos  9'  - y'^,  SIR  e'  ‘ 

* = «'+£. 


As  we  proceed  with  the  summation  of  the  time-  , range-  , and  drift- 
increments  along  the  sound  rays  originally  the  paths  of  their  respective  sound 
wave  surface,  i»e'  approach  the  sound  source,  and  therefore  the  differences  in 
the  X-  and  y~  coordinates  for  points  on  each  ray  at  fixed  levels  approach  a 
minimum.  Hathematically,  this  condition  may  be  expressed  as  follows: 

lim  I-  x^  StR  0 + y^  cos  0 - P + x'^,  stR  0»  - y'^,  cos  0'l<e 

lim  I-  x^  cos  0 - SIR  0 - # + x'^i  cos  0'  e y\i  *•"  0'|<F 

hu  Is  - x'  - £|  < P , 

where  e is  a small  positive  quantity. 

Theoretically  speaking,  the  two  sound  rays  intersect,  since  the  sound 
waves  emanate  originally  from  a point  source.  Practically,  however,  in  the 
process  of  tracing  the  ray  through  the  atmosphere,  utilizing  for  this  purpose 
instrumental  equipment  subject  to  inherent  errors,  it  would  be  a fortunate 
coincidence  were  the  computed  rays  to  intersect.  Therefore  we  resolve  the 
solution  for  the  sound  source  into  the  determination  of  (a)  the  anntmua  dis- 
tance between  two  sound  rays,  and  (b)  the  position  and  orientation  of  the  line 
segment  representing  this  distance  with  regard  to  the  standard  coordinate 
system.  The  sound  source  can  then  be  said  to  be  located  at  some  point  of  this 

line  segment.  The  length  of  the  line  segment  is  an  index  of  reliability  of 

the  sound  ranging  procedure  involved  in  each  instance  of  ranging  operations. 

22.  Minimum  Distance  between  Two  Sound  Rays.— For  the  purpose  of  ready 
identification,  let  the  coordinates  of  general  points  referred  to  the  x,  y,  t, 
system  be  x^  , y^  , on  the  sound  ray  appropriate  to  apace  array  (1)  and 
r , y^  , X,  on  the  ray  appropriate  to  apace  array  (2).  Further,  let  the 
additional  subscript  1 or  2,  be  appended  to  the  coordinates  of  a particular 
point  on  each  ray  to  indicate  the  level  (1)  pr  (2),  respectively,  at  which  the 
point  on  the  ray  is  considered.  Thus,  x^,  , y^^  , x^^  ; x^,  , y^,  , x^,  are 


(!>•) 


(IIT) 


(HR) 


(Itt) 
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th*  coordinatM  of  points  on  ray  (1)  st  the  top  and  bottom  of  s layer  of 


phere. 


Por  the  layer  in  which  t in  equations  (•••)  above  is  a minimum  (at  some 
point  of  which,  therefore,  the  sound  source  is  located),  substitute  for  the 
sound  rays  (1)  and  (8)  a straight  line  Joining  the  end  points  of  the  ray  at  the 
boundaries  of  the  layer.  Provided  the  layer  i*'  not  excessive  in  vertical  dimen- 
sion, this  approximation  is  justified,  since  witAin  a singlt  skallov  laytr  the 
sound  ray  does  not  depart  appreciably  from  a straight  line'\  ro  that  whatever 
errors  are  introduced  are  inconsequential.  It  is  to  be  emphasized,  however, 
that  such  a substitution  Is  not  valid,  in  general,  at  any  part  of  the  sound  ray 
except  the  last  layer  involved  in  the  integration  , since  a slight  error  in- 
troduced elsewhere  along  the  ray  becomes  a cumulative  one  in  the  summation  over 
the  remaining  portion  of  the  assumed  path. 


The  equations  of  the  lines  joining  the  end  points  of  the  rays  in  the 
layer,  expressed  in  two-point  form,  are 


* - *11  _ 

y - = 

* - *11 

< 

*i8-*ai 

)'l8-J'll 

*ia-*ii 

* - *81  _ 

y - >81  _ 

* - *81 

• 

*88~*81 

Xa8~Xsi 

*8 8"* a 8 

(ISO) 


The  denominators  are  proportional  to  the  direction  cosines  of  the  lines, 
, for  line  (1),  X,  , u,  , v,  for  line  (2),  where 


X 

1 


Ai 


(iSI) 


where  + (Xis-yn)*  + ' 

with  corresponding  expressions  for  A,  , Ug  » v,  , Hence  equations  (l»o)  shove 
can  be  written  in  the  symmetric  form  of  space  geometry. 


(U2) 


The  angle  t between  the  lines  is  given  by 


cos  "f  = ^*  + Ut  Ug  + Vg  , 


^^Hor«ov*r,  M art  at  llbartr  to  aub-OlvlOe  anjr  larer  of  constant  a and  Into  as  aanjr  laxtrs  as 

ds 

daslrad  or  raqulrad;  Uiartfora  a sufflelantly  snallow  lajrer  can  alnaxa  be  secured. 


58 


r 


iNgiartiM  ikTMifi  tt* 


Mri  m Awltaatiw  to 


Kto 


or 


(l«>) 


• * (Xj  U,  - X,  Ux)*  (u»  V,  - u,  v^)*  + (v^  X,  - V,  X»)* 

If  the  Wo  lines  given  by  equations  (•*»*)  intersect,  then  their  point 
of  intersection  can  be  re^ily  detereined  by  solving  the  equations  siaultaneously, 
Vb  have,  then. 


1-  Iu  = " - ru  and 


y - y 


IX 


so  that 


+ -^(y  - y^^)  and  * = *,»+  -*(y  - y^)  • 


On  elininating  * and  solving  for  y,  and  subsequently  for  z and  *.  we  finally 
obtain 


*Li-TPe  “ *11 


u. 


InJtl. 


Xj  U,  - X,  Xj  Uj  - Aj  u. 


.itJix- 


Aj  U(  - X,  Uj 

U.  ^1 

Xt  “ Xf  U, 


- >ii)  + 
(**i  - *ii>- 


Xi  Ui 


Xi  U,  - A,  u, 


X«  Vj 


Xi  U,  - Aj  u. 


(j'li  " ^11^ 
■(»si " *ti)  h 
■<^'*1  - >11) 


(IS*) 


The  point  determined  by  relations  (I'*')  above  is  the  intersection  of  line  (1) 
with  the  vertical  plane  containing  the  line  (2);  and  hence  the  designated  sub- 
script. For  the  values  x.  y,  given  in  (US)  for  this  point  of  intersection,  we 
find  on  substitution  in  the  equation 


LUlix  = 

the  vertical  coordinate  »lb  on  line  (3)  when  x and  y have  these  values.  Thus 


*t*  ” *ti 


= *ti  + 


IJ'h  + 


JiiJiiL 


A,  u 


Xi  u,  - A,  X^  u,  - A,  Ui”® 


V. 


Xi  U,  - A,  u. 


•(*21  ~ *11)  “ 


» rv...„(y^^  - >11)  - 

(US) 


Xi  U*  - A,  Ux 


(j'ai  “ ^11^  . 


If  the  two  given  lines  intersect,  then  - vpz  = *ls  • poini 


*Ll  ~ VPS  ' ^L1  -TPS 
sound  source. 


'Ll  - TPS 


determined  above  is  the  position  of  the 


If  the  two  lines  do  not  intersect,  we  can  determine  the  shortest 
distance  between  then.  This  distance  will  be  a line  segment  perpendicular  to 
each  of  the  lines.  The  direction  cosines  A,  v of  the  line  segment  which 
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rcpPMcnto  th«  •hortMt  distwiet  bfttwen  tb*  liaw,  cm  be  }*^ 

of  tbo  dlrMtloo  cooltw*  X , u»  . v»  , Md  x , o,  , v,  . mu*  rroo  m 

h»v«  * 


^ + li,  u + V,  V • 0 


(IM) 


lUainotlng  x >ra  obUin 

u(x,u»  - XjO,)  + v(vjX,  - v,Xj)  = 0 : 
Md  eliolnatlng  u 

x(XjU,  - X,Ut)  + v(u,Vj  - UjV,)  = 0 : 

to  that 


■ A = U = ^ . (1*7) 

lijV,  - u,Vj  VjX,  - v,x»  X»w,  - X.lij 

Substituting  in(l«*)we  have 

t»«*f  = x^^®*''*  " * x*^®*''*  “ 

= V*  t X*  -t  U^(u^v,  - u,vj*  , 
so  that  ve  can  further  vrite 

1 - .1  - X.  = ± 1 , (I**) 

- W,Vi  ■ - v.Xi  ■ X»U,  - X,U»  sin  W 

where  V is  the  angle  between  the  given  lines. 

The  length  d of  the  line  segaent  which  is  the  shortest  distance  between 
the  given  lineb^is  equal  to  the  projection  on  it  of  any  other  segaent  Joining 
a point  orr  one  of  the  given  lines  with  any  point  on  the  other.  Using  points 

*11*  7m  *11  “'I  *ti*  7ai«  *ti*  *'®  ^®7e 

d = X(*ii  - *,4>  + 4(7ii  - 7,i)  + - *,»)  . CW) 


By  use  of  relations  O'**)  we  can  write, 

d = ^ " l^*'^i^^*ii  ” *11^  ^7jXg  ~ ''tXi^(7jj  ~ yfi) 

9%K  T 

■*•  (XiU,  - x,Ui)(»»4  - *,i)l  . 


MrpMdleulw  to  botb. 


60 


A IlMnr  af 


w Iwlloatlw  to 


•Mitai 


or  in  ioteminMt  foni. 


stn  ¥ 


i*xx  ” 

<»ii  " »tt) 

<*»i  - *#i)  ''i 


The  length  of  the  segnent,  then,  is  s function  of  but  only  the  endpoints  of  the 
given  lines  st  the  boundsries  of  the  layer. 


23.  Locus  of  the  Sound  Source. ^Solving  (I**)  simultaneously  uith 
equations  (t*6)  for  a.  u,  v,  -re  have  , 


X = j ° ^ ;u  = ^ ^ ^ ;v  = ^ K “* 


K - ^ 

& = X. 


’ll  « -'ll  'ai  *11  til 


Substituting  from  (ISO)  tre  have  the  simpler  expressions 


sin  f' 


U = ± 


A,  v«  - X,v, 


X.U*  ~ X.U. 


sin  O' 


sin  O' 


where  the  proper  sign  is  chosen  so  that  d is  measured  in  positive  sense  from 

the  given  line  (1)  to  line  (2).  Comparison  of  «li'vp2 

(lOS)  above  will  indicate  the  direction  of  the  perpendicular  required. 


Having  so  determined  d,  the  minimum  distance  between  the  two  lines, 
and  the  direction  cosines  u,  v of  this  common  perpendicular,  we  can  regard 
this  line  segment  as  the  locus  of  the  sound  source.  However,  we  further  re- 
quire the  space  coordinates  of  the  end-points  of  the  common  perpendicular,  so 
that  the  locus  of  the  sound  source  may  be  definitely  known  with  respect  to  the 
origin  of  the  reference  coordinate  system. 


To  determine  the  coordinates  of  the  point  of  intersection  of  the  common 
perpendicular  with  line  (1),  translate  line  (2)  parallel  to  Itself  the  distance  d 
in  the  plane  containing  both  line  12)  and  the  common  perpendicular.  The  x,,  y^, 

X coordinates  of  fixed  points  on  line  (2)  then  become  x^  - \d,  y,  - ud,  x,  - vd, 
and  the  equation  of  line  (2)  in  its  position  after  the  specified  translation  is 


- (x..  - Xd;  _ y - (y,.  - ud)  X - (x 


n 


j 


T I L L E R V 


Solving  nquation  0*»  sinultaMOUslx  vith  the  equation  for  line  (1)  given  in 
(111)  , 

* * - X<i)  + ^ y - ^(y,i  - urf) 


so  that 
or 

>0*0  = + 


= *»i  + - y«)  . 


(^  - ^)  y = . 

u,  * 

U.u. 


('ll  - W)  + J*(>tl  - ui)  - ^ 


— " * Tu  - X U l('.i->>'‘>-'iil%  ^ a 


[(y,i  - u4)  - y»J 


where  *p=o*  *t=o’  *0=0'  '-he  coordinates  of  the  intersection  of  the  perpen- 
dicular wltrh  line  (IK 


given  by 


Lastly,  the  intersection  of  the  common  perpendicular  with  line  (S)  is 


p=p  “ 

*p=o 

+ 

Xd 

II 

II 

y p=o 

+ 

ud 

P=P  “ 

*p=o 

+ 

vd 

2H>  Expressions  of  Position  of  Sound  Source  for  Purpose  of  Msp-Piotting.— 
The  mathematical  solution  for  the  sound  source  is,  therefore, 

1.  Equation  (ISS)  if  the  sound  rays  intersect,  using  as  a criterion 
relation  (ISS)  in  connection  with  (l««)  ; or 

The  common  perpendicular  d,  representing  the  shortest  distance  be- 
tween the  rays,  positioned  as  follows: 

a.  Direction  cosines  u,  v by  equation  (IS2) 

6 Length  d by  equation  (>S0) 

c.  End-points  as  given  in  equations  and  (IBB)  . 

Ve  shall  choose  the  mid-point  of  the  line  segment  d as  the  point 
source  of  sound,  and  utilize  the  magnitude  of  d as  an  index  of 
roliability. 

The  evaluation  of  the  position  of  the  sound  source  in  terms  of  space 
coordinates  permits  the  immediate  adaptation  of  this  determination  to  expression 


62 


m Awl  teat  i«M  to 


A IkMrf  tai*  rmtirtlw  TVaagb  tto 


of  tho  loootloo  of  tbo  oooroo  la  torao  of  /rt4>co*rf<««t«s  for  porpooo  of 
aap-plettla|,  provldod  of  oooroo,  that  tho  grld-ooordiaotoo  of  tbo  orlfta  of 
tbo  roforoaoo  ooordlaato  oyatoa  oro  kaown.  Iho  olovatloe  of  tbo  oooreo  roloi- 
tl«o  to  tbo  orlfla  of  tho  roforoaeo  orotoa  la  dotorainad  io  addltloa. 


Iho  pooltloa  of  tho  aooad  oooreo  aay  alao  bo  oxprooaod  aa  (a)  aa 
oaiolar  doflMtloa  (tra«  •«<■«(*)  of  tho  oooroo  froa  a fliod  poaltloa,  aad 
(O  tho  diataaoo  aloag  tbio  aaiaotb  diroetioa  froa  the  fixed  poaltleh  to  aooroe 
(ee«aeM-reer«).bittiag  C repreaeat  the  true  aiiauth  angle  froa  north,  t the 
range  along  tho  aaiaoth,  m have 


C4. 

= arc  tmi  — 


(IM) 


(X,-  /o)*  + (y,  - fo)*  • 

ehero  (X,,  f,)  aad  (X,i  f,)  are  the  grld-coordiaatea  of  the  sound  souree  aad 
fixed  poaition,  roapootively. 

25.  Solation  for  Norixontal  fropagation  of  Soendy— The  solutions  given  above 
are  applieable  ehanever  the  paths  of  the  sound  waves  considered  are  space 
carves.  The  case  where  the  sound  propagation  is  along  the  horixontal  is  less 
suited  to  rigorous  analysis,  since  irregularities  of  terrain  and  local  nicro> 
aeteorological  influences  both  operate  in  the  direction  of  non-hooogeneity  of 
the  oediun  about  the  sound  path.  Diffraction  over  local  obstacles  adds  to 
conplexity  in  the  interpretatim  of  the  net  effect  of  the  influences  operating. 
Perhaps  the  beat  solation  in  this  case  is  to  assaae  the  net  effect  of  the  dis- 
turbing influences  sero,  and  correct  for  the  observed  azimuth  of  sound  approach 
at  each  array  ^ a factor  proportional  to  the  ratio  of  drift  wind  at  the  ground 
level  to  the  actual  velocity  of  the  sound  along  the  horizontal.  The  intersection 
of  the  two  lines  representing  the  corrected  azimuths  will  then  indicate  the? 
poaition  of  the  sou^  source. 

It  is  to  be  noted,  bomever,  that  if  the  angle  of  inclination  of  the  nave 
noroal  at  a space  array  is  zero,  this  observation  need  not  necessarily  be  in- 
terpreted as  the  horizontal  propagation  of  sound.  It  is  conceivable  that,  on 
occasions,  a particular  sound  wave,^  though  propagated  along  a space  curve  in 
the  ataosphere  and  totally  refract^  at  some  level,  may  arrive  at  the  micro- 
phone array  at  zero  angle  with  respect  to  the  horizontal.  The  sound  wave  aay 
thus  be  considered  to  be  totally  refracted  at  the  kiicrophone  array  by  the 
shallow  layer  of  atmosphere  next  to  the  surface  of  toe  ground.  However,  an  ad»- 
quate  criterion  is  available  in  the  refraction  equation  (37),  since  (a)  if  the 
ground  surface  layer  of  ataosphere  is  the  medium  resultiiig  in  total  refraction 
of  the  sound  wave,  the  inclination  of  the  wave  noraail  at  the  upper  boundary  of 
the  atmospheric  layer  adjacent  to  the  ground  will  be  other  than  n = 1;  whereas 
(A)  if  the  direction  of  sound  propagation  is  alor^  the  horizontal  for  the  en- 
tire path,  then  n = I both  at  the  array  and  at  the  upper  boundary  of  the  sur- 
face layer  of  ataosphere. 

The  solution  for  the  sound  souree  for  case  (a)  is  no  different  from 
that  for  the  general  case  above.  The  solution  for  case  (6)  is,  on  the  basis 
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of  tht  OMuaptlons  obove,  the  angular  correction  e to  be  applied  to  9 where 

“ V. 

e = are  ten  » , 


being  the  velocity  of  sound  along  the  wave  front  nornal.  The  true  azi«ith 
of  the  sound  source  from  the  fixed  position  is  6 + e . 

Using  the  coordinates  indicated  for  the  two  arrays  illustrated  in 
Figure  5t  we  have  for  the  case  of  horizontal  sound  propagation, 

X = y tan  (©  + e)  for  ray  (1) 

S-i  = ton  (e'  e')  for  ray  (2) 

y ' 

But  X = F + X'  , y = K y', 

to  that  X = y ton  (e  + e) 


so  that  finally 


ton  (0'  + e'). 

y - H 


y . p - M ton  -h  p') 

• ton  (0  + c)  - ton  (0'  + e') 

*.  = P - M tan  (0'  -I-  f ') 

ton  (0  + e)  - t»"  (6'  + e') 


ton  (0  + e)  . 


where  x,  , y,  are  the  coordinates  of  the  source  referred  to  the  x,  y,  x, 
coordinate  system. 

Knowing  Fg,the  velocity  of  sound  along  the  normal  and  therefore  in 
this  case  the  same  for  the  horizontal,  and  knowing  the  azimuth  range  S, 
where  R = + y,t  , the  time  of  travel  is  given  by 


t = 


= B. 


If  the  sound  propagation  from  a given  source  to  one  array  is  along 
the  horizontal  and  along  a space  curve  to  another  array, the  solution  for  the 
source  is  determined  as  in  the  general  case  above. 


(1*7) 


(IB«) 


(169) 


(160) 
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In  the  previous  sections  equations  have  been  derived  ubicb  can  be 
utilised  in  tracing  the  sound  ray  through  apace  froa  the  aicrophone  array  back 
to  the  sound  source.  Fbr  this  coaputation  a detailed  knoeledge  of  the  aeteor- 
ological  structure  of  the  ataosphere  is  required,  in  particular,  the  variation 
of  teaperature,  huaidity,  wind  direction,  and  wind  velocity  with  elevation. 

Equations  (s) , (S) , and  (•)  of  Section  II  afford  a aeans  of  deteraining 
the  direction  of  sound  arrival  at  the  aicrophone  array,  both  in  asiauth  and  in 
inclination  to  the  horizontal. 

( 

Equations  (!•)  and  (17)  of  Section  III  express  the  velocity  of  sound  as 
a function  of  teaperature  and  huaidity.  By  this  means  the  variation  of  sound 
veloci^  with  height  can  be  determined  from  a measureaent  of  these  meteorolog- 
ical elements. 

Equation  (37)  of  Section  V is  the  refraction  equation  which  relates  the 
inclination  of  the  sound  wave  normal  at  one  point  with  its  inclination  at  another, 
provided  the  meteorological  elements  at  each  point  are  known.  ^ means  of  this 
equation  it  is  possible  to  trace  the  sound  ray  through  the  ataosphere  to  its  point 
of  total  refraction.  Equation  (ss)  expresses  the  time  of  travel  of  the  sound  ray 
through  a layer  of  atsK>sphere,  and  equations  (3s)  and  (32)  or  (35)  afford  the  means 
of  computing  the  Ineteorological  corrections"  along  the  path  of  propagation. 

Equations  (tss),  (IS5),  (ISO),  (152),  (I5s),  and  (l55)  of  Section  VII 
afford  a means  of  computing  the  position  of  the  sound  source  matheaiatically 
from  the  space  curves  representing  sound  paths  to  two  microphone  arrays. 

In  a subsequent  report  a description  of  a practical  field  procedure  for 
rapid  calculation  of  the  sound  path  will  be  presented.  A special  slide  rule, 
with  scales  incorporating  the  features  of  the  above  equations,  will  be  described, 
and  its  operation  for  compulation  of  meteorological  corrections  for  sound  ranging 
will  be  discussed. 
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